THE ENGINEER GRAPPLES WITH NONLINEAR PROBLEMS! 
THEODORE VON KARMAN 


I do not believe that one could connect justly the name of Gibbs 
with practical applications of applied mathematics, for his main in- 
terest was certainly centered on basic conceptions of mathematical 
physics. Nevertheless, for example, his beautiful work on graphical 
methods in thermodynamics is a brilliant example of the presentation 
of theoretical relations in a form which appeals to the engineer. 

This lecture is intended as an effort to improve the convergence be- 
tween the viewpoints of mathematics and engineering. Thus, I feel 
it is not inappropriate to dedicate it to the memory of Josiah Willard 
Gibbs. 

Engineering mathematics is generally considered as a collection of 
mathematical methods adapted for the solution of relatively simple 
problems. These problems often might require lengthy numerical cal- 
culations or graphical constructions, but supposedly can be worked 
out without the use of advanced methods of mathematical analysis. 
This description was perhaps correct some decades ago; today a large 
group of scientific workers is engaged in applying various methods of 
classical and modern analysis to problems in electrical, civil, mechani- 


cal, aeronautical and also chemical engineering. It is not possible to 
give an exhaustive list of all types of problems which require the ap- 
plications. of advanced analytical methods. In the following table 
merely some of the most important engineering problems and the 
mathematical concepts and methods involved in their treatment are 


indicated: 


ASSOCIATED TOPICS OF ENGINEERING AND MATHEMATICS 


Mathematical topics Engineering problems 


Vector algebra, systems of linear Engine dynamics, vibration of sys- 


equations. Tensors and matrices. Alge- 
braic equations. Ordinary differential 
equations with given initial conditions. 
Elementary operational calculus. 


tems with a finite number of degrees of 
freedom. Rotating electric machinery. 


Ordinary differential equations and 
their boundary problems. Eigenvalues 
and eigenfunctions. Expansion in or- 


Equilibrium, buckling and harmonic 
vibrations of beams. Critical frequencies 
and speeds. One-dimensional problems 


1 The fifteenth Josiah Willard Gibbs Lecture, delivered at Columbus, Ohio, De- 
cember 27, 1939, under the auspices of the American Mathematical Society with the 
cooperation of the American Association for the Advancement of Science. 
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Mathematical topics 


Engineering problems 


thogonal functions. Calculus of varia- 
tions of simple integrals. 


of stationary flow of matter, ele tricity 
and heat. 


Functions of complex variables. Con- 
formal transformation. Logarithmic po- 
tential. 


Two-dimensional velocity, electro- 
magnetic and temperature fields. 


Vector and tensor analysis. Partial 
differential equations of the elliptic type 
and their boundary problems. Expansion 
in orthogonal functions. Integral equa- 
tions. Calculus of variations of multiple 
integrals. 


Stress distribution in elastic bodies; 
equilibrium, buckling and harmonic vi- 
bration of plates and shells. Three-di- 
mensional flow and field problems. 


Partial differential equations of the 
parabolic and hyperbolic type. Advanced 
operational calculus. 


Transient phenomena of heat con- 
duction and wave propagation (variable 
flow in channels, water hammer in pipes, 
surges in cables, waves in solids and 
fluids, electromagnetic waves). 


Though most of the various branches of mathematical analysis ap- 
pear in this table, it does not follow that the engineer has to master, 
for example, the complete theory of linear differential and integral 
equations or the calculus of variations. He is mostly interested in 
methods of obtaining approximate solutions. However, it is interest- 
ing to note that the methods which are perhaps the most powerful 
tools of the mathematician for proofs of existence of solutions have 
been proved to be the most helpful in obtaining numerical approxima- 
tions, for example, step-by-step integration, iteration, successive ap- 
proximation, direct methods of finding maxima and minima. The 
iteration method is known to most engineers as the Stodola or 
Vianello method. The direct methods of the calculus of variations 
are known as the energy method, Rayleigh-Ritz method, Galerkin 


method, and so on. 


The problems listed in the above table provide a happy hunting 
ground to the mathematician who is interested in applying his science 
for practical purposes. We know that not all mathematicians have 
this interest; a great British mathematician is credited with the re- 
mark that his field—the theory of numbers—is the only field of 
mathematics that, fortunately, has not been prostituted by being 
useful to some practical end. Others believe that the engineering ap- 


plications of mathematics do not require ingenuity in the mathemati- 
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cal sense; however, there are many engineering problems that are far 
from being an easy mark for the hunter. Most of the problems com- 
piled above are readily accessible to the methods of mathematical 
analysis because they are linear problems. In many cases the engineer 
can linearize his problem by means of simplifying assumptions and a 
mathematical text will easily supply him with all the help he needs. 
However, if the engineer has a real nonlinear problem, that is, one 
which loses its sense by linearization, very often he has to grapple 
with it by himself. 

In this paper the attempt is made to give a review of nonlinear 
problems which are of practical interest and for whose solution meth- 
ods have been worked out or at least suggested by mathematicians, 
physicists or engineers. The review is not exhaustive, but it was at- 
tempted to include examples of the most typical problems. The ques- 
tion of existence and uniqueness of solutions is only superficially 
+yuched; emphasis is laid on methods for obtaining approximate solu- 
tions. Many physicists and engineers have the feeling that so-called 
mathematical proofs of existence and uniqueness of solutions very 
often merely verify facts which are more or less evident from physical 
considerations. Unfortunately, in most nonlinear problems physical 
reasoning is not sufficient or not fully convincing, so that in these 
cases the questions of existence and uniqueness represent a real chal- 
lenge to the mathematician. Certain investigations of H. Poincaré 
and T. Levi-Civita can be mentioned as examples in which math- 
ematical reasoning has virtually extended the range of physical vi- 
sion. The modern methods of function spaces are likely to make 
essential contributions to the problem of existence and uniqueness 
of solutions of nonlinear equations and their boundary problems. 
Promising steps in this direction have already been made. However, 
the writer is afraid that this aspect of the subject is beyond his own 
competency and the scope of this paper. 


1. Nonlinear oscillations. General remarks. The theory of mechan- 
ical oscillations commonly used in engineering practice is based on 
the assumption that the restoring forces are proportional to the de- 
flections and the damping forces to the velocities, that is, to the first 
derivatives of the deflections with respect to time. If we assume con- 
stant masses, the equations of motion of an oscillating system are 
linear, containing as coefficients constant inertia, damping and spring 
factors. We obtain analogous equations for the oscillations of an elec- 
tric circuit by assuming constant impedances, that is, by assuming 
that the network contains only constant reactances, resistances and 
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capacities. The oscillations governed by linear equations may be 
called “linear oscillations.” 

Let us consider the oscillations of a single mass m assuming that 
the deflection x produces a restoring force —kx and the damping 
force is equal to —8z, where t =dx/dt=v is the velocity of the mass. 
Then the equation of motion for the free oscillations of this mass is 


(1.1) mi + kx + Bx = 0. 


Introducing the natural frequency w of the undamped oscillations, 
that is, the frequency for the case 8B =0, we have 


(1.2) 
m 


We will generalize equation (1.2) in the form 
(1.3) + wx = f(x, 


This general form of the equation of free oscillation includes the 
cases of nonlinear damping and nonlinear restoring force. Often the 
deviation from linearity causes only a small deviation of the linear 
oscillation without changing the general character of the motion. For 
example, the simple harmonic oscillation, which we would obtain in 
the linear case, will be accompanied by harmonics of small amplitude 
if the deviation from linearity is taken into account. However, in 
some cases the whole character of the oscillatory motion changes. 
Slight negative damping in the range of small deflections may cause 
“self-excited periodic oscillations,” whose period is quite different 
from the period of the undamped harmonic oscillation. Nonlinearity 
of the spring factor may fundamentally change the resonance phe- 
nomena, that is, the response of the system to external periodic forces. 
In the case of a pure linear system the response to a periodic force 
consists of a forced oscillation whose period is equal to that of the 
force and of free oscillations, whose period is the natural period of the 
system. If the two frequencies coincide, resonance occurs. One can 
easily visualize that in a nonlinear case the periodic force will excite 
also higher harmonics, that is, oscillations whose frequencies are mul- 
tiples of its own frequency. It is more difficult to see offhand that in 
such a system an external force can also excite oscillations whose 
frequency is a fraction of its own frequency. This phenomenon is 
known as “subharmonic resonance.” The mathematical analysis of 
the phenomena of self-excited oscillations and subharmonic resonance 
makes use of certain methods developed by H. Poincaré in connection 
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with his astronomical investigations and shall be discussed briefly in 
the next two sections. 


2. Self-excited nonlinear oscillations (relaxation oscillations). Let 
us assume that the function f(x, <) in equation (1.3) has the form 
f(x, %) =ub(x)z, where yu is a positive constant parameter. It is seen 
that in this case the damping is proportional to the velocity v=, but 
the magnitude of the damping factor is in general a function of the 
deflection. The equation (1.3) in this case takes the form 


(2.1) + = ph(x)x. 


Equation (2.1) can be reduced to a differential equation of the first 
order by considering =v as function of x. We obtain 


dv x 
(2.2) — = — w— + po(x). 
dx v 


This equation can be solved by the method of isoclines. Consider first 
the case 1 =0. The integral curves of (2.2) are in this case concentric 
ellipses, given by 


(2.3) v? + wx? = const. 


One of these ellipses is shown in Fig. 1a. Since dt =dx/v, we obtain 
the x, v values corresponding to the actual motion by proceeding 
clockwise along the integral curves. The corresponding motion is a 
harmonic oscillation with constant amplitude. 

Let us now assume that $(x) =1, that is, the damping factor is con- 
stant. It is known that the resulting motion is a harmonic oscillation 
with decreasing amplitude if «<0 and increasing amplitude if »>0. 
We obtain logarithmic spirals, which in the stable case (u<0) con- 
verge toward the origin (Fig. 1b) and in the unstable case (u>0) di- 
verge from the origin (Fig. 1c). 

We shall now assume that (x) is variable and investigate the exist- 
ence of periodic solutions. In the x, v plane a periodic solution ap- 
pears as a closed curve. Let us multiply both sides of equation (2.2) 
by vdx and integrate along the closed curve. Then we see that 
US o(x)vdx =n $ o(x)v*dt =0. Hence ¢(x) must change its sign during 
the motion. In other words, to obtain a periodic motion it is a neces- 
sary condition that the damping factor change its sign as function of 
the deflection. 

Van der Pol studied in detail an interesting special case assuming 
that ¢(x) =1—x?. It is seen that for x <1 we have negative damping 
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which becomes positive when x>1. Fig. 2 and Fig. 3 show integral 
curves in the x, v plane for 4.=0.1 and w=1.0 respectively. The in- 


Fic. 1a, b, c 


tegral curves near the origin have the shape of logarithmic spirals, 
but instead of diverging to infinity they converge to a closed curve. 
The integral curves originating at infinity converge to the same closed 
curve, which represents a periodic motion with constant amplitude. 
For small values of u, for example, 4=0.1, the closed curve is only 
slightly different from an ellipse, but its shape varies quite radically 
with increasing yu. The corresponding motions for »=0.1, 1.0 and 10 
are shown in Fig. 4 where the deflection is plotted as function of time. 
It is seen that for »=0.1 the motion is very nearly harmonic, whereas 
for 4=10 the oscillation is made up of sudden transitions between 
deflections of opposite sign. The period of the harmonic motion is de- 


v y 
(a) 
\ | 
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termined by the spring factor and the mass. The period of the motion 
corresponding to 4 =10 is much larger and depends both on w and yp. 


Fic. 2 


From B. Van der Pol, Philosophical 
Magazine, vol. 2 (1926) 


For this reason Van der Pol calls this type of motion a relaxation 
oscillation. The term self-excited oscillation is also often used. The 
physically interesting point is that we obtain a periodically undamped 
oscillation in a damped mechanical or electrical system because the 
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damping coefficient is variable and for small deflections has a slight 
negative value. 
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From B. Van der Pol, Philosophical 
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The stability of this type of motion can be investigated by the 
method of the “equations of variations.” Let us assume that 


(2.4) (x, Xo, 0%) 
x 
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From B. Van der Pol, Philosophical 
Magazine, vol. 2 (1926) 


is a solution of (2.2) passing through the point xo, v9. Then we have 


(=) == d\\~ @ 


or 
d av wx dv 
dx \0xo v2 Oxo 


where v is substituted from equation (2.4) as function of x, xo and 2p. 
It is seen that by integrating between the limits x; and x2 we have 


dv 2 wx 1 dv #2 
[8 =| f —Fartuf dx. 


1 1 1 


Applying this equation to a closed cycle, we obtain 
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In a similar way we obtain 


or 


Ov ov 


Let us consider now an integral curve in the neighborhood of the 
closed cycle passing through the point xo, vo+ 462, and denote the 
value of v obtained after completing one turn along the closed cycle 
by vo+ 

Then according to (2.9) 


(2.10) = 50; exp {ug 


It is seen that if $o(x)dt>0 the deviation of a neighboring integral 
curve from the closed cycle increases; if $@(x)dt<0, decreases for 
each turn. We conclude that the motion corresponding to a periodic 
solution x(t) is stable if $(x)dt <0. 

The most important applications of the theory of oscillations of 
systems with nonlinear damping sketched in this section are in the 
field of radio. Both the theory and its applications have been treated 
by several authors. An extensive list of references is given in the paper 
of Van der Pol mentioned in the Bibliography at the end of this papér 
[5]. Van der Pol pointed out that in addition to oscillations in me- 
chanical and electrical systems many other periodic phenomena have 
the character of relaxation oscillations. He found that differential 
equations similar to those of the relaxation oscillation describe with 
fair approximation various periodic biological processes, for example, 
the “struggle for life” of certain animal groups investigated mathe- 
matically by A. J. Lotka and Vito Volterra. Van der Pol believes that 
even periodic business cycles show a certain analogy to the relaxa- 
tion oscillation of a physical system. The essential condition for such 
oscillations is negative damping for small deviations and a rather 
rapidly increasing positive damping for large deviations from the 
equilibrium position. The psychological response of certain groups of 
people to changing business conditions shows doubtless some anal- 
ogy to the behaviour of mechanical systems capable of relaxation 
oscillations. 


| 
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3. Subharmonic resonance. Let us now assume that the function 
f(x, 2) in equation (1.3) has the form f(x, z) = —y[z+y(x) ]; that is, 
the system has linear damping but the elastic restoring force is a 
nonlinear function of the deflection. The center of interest is the re- 
sponse of the system to the action of a periodic force; as we men- 
tioned, it was observed that in many cases the oscillations of the 
natural frequency w of the system are excited by a periodic force 
whose frequency is a multiple of w. This curious phenomenon is called 
“subharmonic resonance.” 

If we consider first the linear system corresponding to 4.=0 and 
assume that a periodic force a sin nwt is applied to the system, the 
general solution of the equation of motion is 


(3.1) x c® sin (wt — a) + ol a) sin mot . 

The first term on the right side represents a free oscillation of the 
system with an arbitrary amplitude c\ and phase angle a, and 
the second term a forced oscillation whose frequency is equal to the 
frequency of the impressed force. If we introduce slight constant 
damping [u=const., ¥(x)=0] the free oscillation dies out and the 
forced vibration remains with a slightly modified amplitude. Let us 
now assume that ¥(x) #0. Several authors reached the conclusion 
that in this case the solution consists of higher harmonics with the 
frequencies 2nw, 3nw, - - - . In fact, if one would start from the forced 


oscillation 
a 


w?(i — 


sin 


and determine a next approximation by substituting # in (x), it 
is easily seen that one would obtain forced oscillations whose fre- 
quencies are multiples of mw. However, there is a fallacy in this con- 
clusion because in order to obtain the most general type of periodic 
motion we have to start with the complete solution (3.1) and put 


(3.2) x= ye) 4 , 


The function x has two undetermined constants c® and a°. They 
must be chosen in such a way that x is a periodic function of ¢ with 
the period T=27/w. Substituting (3.2) in the differential equation 


(3.3) + wx = pf(x, + asin not 


and equating terms of the order u on both sides of the equation we 
obtain 
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(3.4) ED + = f(x, 


The solution of (3.4) that satisfies the initial conditions x =z =0 
for t=0 is given by 


t 
(3.5) = f f[x(r), £(r)] sin w(t — r)dr. 
0 
This solution corresponds to a periodic motion with the period 27/w, if 


= — x(x), (r)] sin wrdr = 0, 
(3.6) rf 
-f f[x(r), cos wrdr = 0. 


Equations (3.6) are two equations for the undetermined parame- 
ters c® and a, Since x(t) and #(t) are periodic functions of ¢ 
with the period 27/w, we can expand f[x‘(#), #(t)] in a Fourier 
series of the form }>,a, sin nwt +)°.b, cos nwt. According to equations 
(3.6) the coefficients a; and b; vanish. Therefore, the solution of equa- 
tion (3.4) will have the form 


2 sin nwt bo 
— n?) w? 
(3.7) 
= Cos Nwt 


+> 


n=2 w(1 


Substituting x=x©-+ypx™ and into f(x, in equa- 
tion (3.3), we obtain an equation for x and two equations analogous 
to equations (3.6). If we neglect higher powers of u, we have two lin- 
ear equations for c and a. Thus we are able to determine x as 
function of ¢ to any desired accuracy by successive approximations. 

It is seen that the final solution will contain a harmonic oscillation 
with the frequency w. This oscillation is called the subharmonic reso- 
nance. It has a definite amplitude c and phase a, which can be de- 
termined by the method of successive approximations described 
above. Subharmonic resonance occurs if the equations (3.6) have real 
roots c and a. Also it is necessary that the determinant of the 
linear equations for c™ and a” be different from zero. 

The physical meaning of the conditions (3.6) is the following. The 
expression pf(x, <) is an approximate expression for the force re- 
sulting from damping and elastic restraint. If the conditions (3.6) 
are satisfied, that is, the harmonic components of the force with the 
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frequency w are zero, the work done by this force on the deflection 
c® sin (wt —a™) vanishes. 

The problem of subharmonic resonance was treated by a great 
number of authors. The most important papers are mentioned in the 
bibliography appended at the end of this lecture. Tuckerman also 
constructed a simple and ingenious model for the demonstration of 
subharmonic resonance.” His attention was called to the problem by 
the fact that certain parts of an airplane can be excited to violent 
oscillation by an engine running with a number of revolutions 
much larger than the natural frequency of the oscillating parts. 

The method of this section can be applied to forced oscillations of 
systems with nonlinear damping treated in §2. It has been found that 
also in such systems a force of the frequency mw can produce harmonic 
oscillations of the frequency w. This procedure in the radio-technique 
is called “frequency demultiplication.” Frequency demultiplication 
can occur in systems with positive or negative damping for small de- 
flections, that is, also in such systems which are not capable of self- 
excited relaxation oscillations. 

The mathematical problems outlined in §§2 and 3 represent simple 
cases of problems treated by H. Poincaré in his various researches. 
The mathematician will readily recognize in the closed curves shown 
in Figs. 2 and 3 the “cycles limites” studied by Poincaré in his paper, 
“Sur les courbes définies par une équation diffeientielle,” Oeuvres, vol. 
1. The method used in this section is analogous to the methods of the 
perturbation theory given in the Méthodes nouvelles de la mécanique 
céleste. = 


4. Nonlinear problems in the théory of elasticity. Large deflections. 
In the theory of elasticity usually two fundamental assumptions are 
made: small deflections and linear stress-strain relations (Hooke’s 
law). Under these assumptions the problem of stress distribution in 
elastic systems is governed by linear partial differential equations. 
However, if we give up one or both of these assumptions, we obtain 
nonlinear equations. We consider in this section the case in which 
Hooke’s law is assumed to be correct, that is, the stresses are linear 
functions of the strain components, but the deflections produced by 
external loads are so large that in the expressions for the strain com- 
ponents the squares and products of the deflections and their deriva- 
tives cannot be neglected. A classical example of this type is the 
problem of the “Elastica,” that is, the bending and twisting of thin 


2 The author is indebted to Dr. Tuckerman for an extensive list of references and 
for a film shown during the presentation of this lecture. 
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rods. Many of the great mathematicians of the eighteenth and nine- 
teenth century, for example, James and Daniel Bernoulli, Euler, 
Lagrange and Kirchhoff, were interested in this problem. 

Let us consider a simple case: two equal and opposite forces of 
* magnitude P, acting along the x axis, are applied to the endpoints of 
a uniform strut. The endpoints can move freely along the x axis, but 
are restrained to this axis. The strut can be deflected in the x, y 


P/P; 


2.0 


1.0 


0 025 0.50 
Wmax/ I 


Fic. 5 


plane. Then we obtain for the deflection w in the y direction the non- 
linear differential equation 


= 
ds? B ds 
where s is the length of the arc measured along the center line of the 
strut and B is its bending stiffness. We assume that a positive value 


of P corresponds to compression in the strut. 
For small deflections equation (4.1) becomes 


(4.2) —+—w= 


If we consider P as parameter in equations (4.1) and (4.2), the 
eigenvalues of (4.2) which comply with the boundary conditions w =0 
for x=0 and x=/ (I the length of the strut) give us the points of 
bifurcation of the equilibrium (that is, the branch points of the solu- 


= 
0. 
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tions). These eigenvalues are equal to P,=Bn?x?/]?, where n is a 
positive integer. The first characteristic value (n=1) determines the 
so-called buckling load, or “Euler-load,” Pg=2?B/I? of the strut. The 
eigenfunction corresponding to Pg gives the shape of the buckling 
column leaving the magnitude of the deflection undetermined. If we 
are interested in the force-deflection relation beyond the buckling 
load, we have to integrate equation (4.1), which can easily be done 
by use of elliptic functions. Fig. 5 shows the ratio between the force P 
and the buckling load Pg as function of the ratio between the deflec- 
tion Wmax Of the midpoint and the length / of the beam. 

The analogous but much more complicated problem of the finite 
deflection of plates loaded beyond their buckling limits has great im- 
portance for the determination of the strength of certain aircraft 
structures. The so-called stressed skin (or monocoque) method of con- 
struction widely employed for metal airplane wings and fuselages is 
based on the idea of loading the thin metal skin beyond its buckling 
limit. In fact the skin in the buckled or “wave” state is able to carry 
stresses which in many cases are large multiples of the buckling stress. 
Then the maximum load is determined by the ultimate strength of 
the buckled skin and the aim of theory is the determination of the 
stress distribution and the maximum stress occurring in the wave 
state. 

‘To fix the idea let us assume that a plate of constant thickness ¢, 
which in the undeflected state is parallel to the x, y plane, is sup- 
ported freely along two parallel edges, y= +5, and is subjected to 
uniformly distributed compressive forces acting normal to its other 
two edges, x= +a. The problem.is much more complicated than the 
problem of the Elastica for the following reason: In the case of the 
deflected thin rod, the resultant of the normal stresses in an arbitrary 
cross section is determined statically by the end loads. In the case of 
the deflected plate, the distribution of the resultant stresses acting 
in the plate depends on the distribution of the normal deflection 
w(x, y). Let us denote the components of the resultant stresses per 
unit length of the plate by az, o, and 7,,. It is known that these stresses 
are in equilibrium if they can be derived from a function F(x, y)— 
called Airy’s function—by the relations 


=—» = 

Ox? 
If the squares and products of the deflections and strain components 


are neglected, the function F(x, y) satisfies the linear differential equa- 
tion 


(4.3) 
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(4.4) AAF = 0. 

On the other hand if no normal external load is applied to the plate, 
the normal deflection w(x, y) obeys the equation 

(4.5) BAAw = 0 

where B is the bending stiffness of the plate per unit length. However, 
if we take into account the quadratic terms in the deflections and 


strain components, we obtain the following two simultaneous equa- 
tions for the functions F(x, y) and w(x, y): 


dw 
Oxdy Ox? dy? 


Oy? Ox? = Ay’? Oxdy Axdy 


(4.6) 


These equations were first given by the present writer in 1910. 
Let us assume that in the undeflected state 


o,= —Il, o,=7,,=0. 


Then substituting these values in the second equation (4.6) we obtain 
the linearized equation for w(x, y): 

(4.7) BAAw + II, — = 0. 

Ox? 
This equation has solutions different from w=0 which satisfy homo- 
geneous boundary conditions for certain eigenvalues of the parame- 
ter II. 

For example, assuming hinged support at the edges x= +a, y= +), 

the eigenvalues will be given by 


B xa? 
(4.8) n == +2) 


where m and m are integers. These eigenvalues represent the branch 
points of the equilibrium for the plate. The lowest eigenvalue gives 
the buckling load IIz. 

If the load II is only slightly higher than this lowest critical value 
IIz, the system of equations (4.6) can be solved by successive ap- 
proximations. However, if the stress applied is many times larger 
than the “Euler stress” IIz, this method converges much too slowly 
for practical purposes. A satisfactory solution of equations (4.6) for 
large values of II/IIz has not as yet been given. The present author 
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in a joint paper with E. E. Sechler and L. H. Donnell attacked the 
problem with a somewhat crude approximation using the notion of 
the “effective width.” The simplifying assumption was made that 
two narrow strips near the supported edges carry essentially the whole 
load. The computations of H. L. Cox, M. Yamamoto and K. Kondo, 
and K. Marguerre are more elaborate; however, they also are not 
devoid of arbitrary assumptions. Marguerre uses the Rayleigh-Ritz 
method for the approximate solution of equations (4.6). Since he re- 
tains only a few terms of a trigonometric series, his solution if II/Iz 
is very large does not approximate very well the shape of the deflected 
plate observed in the experiments. 

Recently K. Friedrichs and J. Stoker attacked the problem of stress 
distribution in a circular plate beyond the first buckling, assuming 
axial symmetry for the deflection. Due to this assumption, the prob- 
lem is simpler than that of the rectangular plate. They were able to 
integrate the system of equations (4.6) which in this case are reduced 
to ordinary differential equations. They investigated the solution for 
the whole range of values of the ratio II/IIz, where II is the radial 
stress per unit area applied on the circumference of the plate and Iz 
is the first critical value of II. Below the buckling load, that is, for 
II <IIz, the circumferential stress o, is uniformly distributed over the 
whole plate and is equal to II. If II exceeds the buckling stress, o, in- 
creases near the edge of the plate and decreases in the center part. 
For large values of II/IIz, o. in the center part becomes negative; it 
changes from compression to tension, whereas the whole external load 
is practically carried by a narrow strip near the edge, which is analo- 
gous to the effective width introduced by the present writer in the 
theory of the rectangular plate mentioned above. In this limiting case 
the method of perturbations can be applied with success since the 
perturbation term contains the derivative of the highest order. This 
amounts to an asymptotic integration of the differential equations re- 
sulting from equations (4.6). Let us denote the average value of a. 
by p and the edge stress by p.. Then Friedrichs and Stoker obtained 
the following asymptotic approximation for the relation between p 


and 


where a is the radius of the plate. 
The same relation with a numerical factor of the same order of 
magnitude is easily obtained by the method which the present author 
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used in the case of rectangular plates. Hence, the analysis of Fried- 
richs and Stoker can be considered as a verification of the authors’ 
simplifying assumptions. 

The theory of large deflections of plates subjected to lateral loads 
requires the integration of the nonhomogeneous system of equations 


Oxdy Ox dy? 


Oy? Ox? dx? dy? Oxdy Axdy 


(4.10) 


] = p(x, 


where p(x, y) is the lateral load per unit area. This system of equa- 
tions can be solved by successive approximations. Denoting the solu- 
tion of the nonhomogeneous linear equation 


(4.11) AAw = p/B 


by wo(x, y), we substitute this function in the first of the equations 
(4.10) and determine a first approximation of Fi(x, y); then, substi- 
tuting the derivatives of F; in the second equation (4.10), we calcu- 
late the next approximation w,;(x, y), and so on. The system of 
equations (4.10) was solved for certain simple contours and loading 
conditions by S. Timoshenko and S. Way. 

The problem of large deflections of curved rods was treated by sev- 
eral authors. Thus R. Grammel investigated the deformation of a cir- 
cular ring by uniformly distributed moments along the circular ring 
axis, which tend to turn the cross sections of the ring around their 
centers. He found that at a critical value of the moment, the ring 
jumps over into a new equilibrium position. 

Another interesting case is the deflection of a flat arch under lateral 
load where the deflection is of the same order of magnitude as the 
height of the arch. Let us assume, for example, that the two endpoints 
of a flat circular arch are pin-joined at two fixed points on the x axis 
and the arch is loaded by a uniform normal load p per unit length. 
If R is the radius of curvature of the arch, the uniform normal load 
produces a compressive force equal to pR and the first bifurcation 
point of the equilibrium is obtained approximately by Euler’s equa- 
tion for a straight strut: 


B 
(4.12) (pR)z = ts 


where B is the bending stiffness, assumed to be constant along the 
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length and / is the length of the arch. The nonlinear equation, which 
takes into account the large deflections, is analogous to the equation 
of the Elastica mentioned before. However, whereas in the case of the 
originally straight Elastica the thrust to be applied increases with the 
deflection, in the case of the originally curved arch, we find that the 
load p which is necessary to keep the rod in the deflected shape can 
be smaller than the critical load corresponding to equation (4.12). 
Thus, as Fig. 5 shows, a straight beam under axial load has only one 
equilibrium position for values of the parameter P, which are smaller 
than the first eigenvalue Pz, whereas the flat arch has several possible 
equilibrium positions for values p< pg as shown in Fig. 6. Hence the 


P/Ps 


05 


A 


— > 
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engineer is more interested in the smallest value of the parameter , 
which is compatible with the existence of such equilibrium positions 
than in the value pz corresponding to the branch point of these solu- 
tions. To be sure a continuous transition from the unloaded position 
to such an equilibrium state is only possible if the load p first reaches 
the critical value pg. However, imperfectness of the structure and 
slight deviations from the uniform load distribution might cause a 
jump between equilibrium states along AB and BC, practically at 
any value p between Pmin and px. 

It seems to the writer that these considerations are of great im- 
portance for the explanation of certain systematic discrepancies 
found between the breakdown loads of curved shells observed ex- 
perimentally and computed analytically from the eigenvalues of line- 
arized equations. For example, the differential equation for the small 
deflections of a thin-walled spherical shell of constant thickness sub- 
jected to uniform external pressure was carefully investigated by sev- 
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eral authors, for example, R. Zoelly, E. Schwerin, and A. Van der 
Neut. The critical pressure corresponding to the smallest eigenvalue 
of the differential equation is found to be equal to 


2 


where E is Young’s modulus, v Poisson’s ratio, ¢ the wall thickness, 
and R the radius of the spherical shell. The experiments show that 
the spherical shell collapses under a pressure that is only one-quarter 
to one-third of this value. Also, the shape of the buckled shell is very 
different from the wave form predicted by the theoretical eigenfunc- 
tion that corresponds to the smallest eigenvalue. In a paper published 
jointly with H. S. Tsien the writer has shown by approximate calcula- 
tion of the large deflections the existence of more than one equilib- 
rium state at values of the parameter p which are much smaller than 
the value given by equation (4.13) and are of the order of the values 
observed experimentally. A similar discrepancy between the predic- 
tion of the linear theory and the experimental results in the case of 
the buckling of cylindrical shells under compression applied parallel 
to the cylinder axis probably can also be explained by analogous con- 
siderations. 

Let us finally consider the problem of large deformations of a three- 
dimensional isotropic elastic continuum. This problem, although 
treated in great detail in the past by several authors, was more or 
less neglected until interest was stimulated by the possible applica- 
tion of tensor calculus as shown by recent publications of E. Trefftz, 
F. D. Murnaghan, and M. A. Biot. The first treatment is due to 
G. Kirchhoff (1852); later J. Boussinesq (1872), M. Brillouin (1891) 
and J. Finger (1894) were interested in the problem. An excellent 
presentation of the subject can be found in a review of E. and F. Cos- 
serat (1896). The analysis of all these authors is based on the funda- 
mental assumption that the state of the matter at a point P, whose 
coordinates before the deformation are x;, is determined by the de- 
formation of the infinitesimal neighborhood of P, that is, by the co- 
efficients of the linear transformation 


0%; 
(4. 14) dk; == d(x; + dx, 
Ox, 
where x; and &; are the coordinates of P before and after deformation. 


The vector whose components are #;—x;=£; is the displacement vec- 
tor. The transformation (4.14) has nine coefficients. However, a rigid 
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rotation does not produce stresses, and therefore the stresses can only 
depend on the coefficients of the square of the line element dz;. We 
write 


(4.15) dk; = gudxidxy 
where 
(4.16) Lik = 


and the a’s are the coefficients of the transformation (4.14). In the 
case of rigid rotation, gi. = 5; (that is, if i=k and if iXk). 
Hence, we conclude that the six quantities g;,— 6; determine the pure 
deformation of an element. E. and F. Cosserat, and later E. Trefftz, 
defined these quantities as the components of the strain €; and con- 
cluded that the stress-strain energy of an isotropic body must be a 
function of the invariants which can be constructed of the quanti- 
ties €;,. If the form of the energy function is assumed, the components 
of the stress tensor are obtained by variation of the energy and ap- 
plication of the principle of virtual displacements. 

If the squares and products of the derivatives of the te S are neg- 
lected, giz = 55 +0£;/Ox,+0&,/Ox;. Hence, in this case the transforma- 
tion (4.14) can be split into two parts: 


1 / 0; 
(4.17) dé; = + =)a Xk +— ~(= = =) dx,. 


The first term represents an .infinitesimal pure deformation, the 
second term an infinitesimal rigid rotation. The quantities €, 
= 4(0§;/dx,+0,/0x;) are known as components of the strain ten- 
sor. If Hooke’s law is valid, the stress-strain energy is a homoge- 
heous quadratic form of the €’s, and the components of the stress 
tensor are derivatives of this quadratic form with respect to the €;;’s. 

If higher terms in the derivatives of the £,;’s are taken into account, 
two difficulties arise: 

(a) The separation of the pure deformation from the rotation and 
the definition of the strain tensor becomes somewhat arbitrary. 

(b) One must differentiate between stresses referred to the surface 
elements of the undeformed and deformed system. 

M. A. Biot suggested the following definition of the rotation and 
the strain tensor: The linear transformation 


(4.18) i= 


Xk 
Xk 
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can be replaced by two successively applied linear transformations, 
where the first one: 
(4.19) xi = 


is an orthogonal transformation representing a finite rigid rotation 
and the coefficients of the second one: 


(4.20) = 


constitute a symmetrical matrix, that is, yuw=‘i. Then the first 
transformation defines in a unique way the rotation and the second 
one the pure deformation. The strain tensor is defined by 


(4.21) = — 


We notice that the yx’s are equal to the gx’s only in the first ap- 
proximation; they differ if higher terms in the derivatives of the &;’s 
are taken into account. 

F. D. Murnaghan’s definition of the strain tensor is based on the 
inverse transformation 


Ox; 


(4.22) dx; = — 


OX; 
If we express the square of the line element ds; by d5;, we obtain an 
expression of the form 


(4. 23) ds; = 


Then we can consider 


(4.24) €ik = hi) 


as the components of the strain tensor. 

Using different conceptions for the strain tensor, Hooke’s general- 
ized law appears in different forms in the publications mentioned. 
Fortunately the question is of academic rather than practical impor- 
tance. With the exception of a few materials with large elastic exten- 
sion, such as rubber, the strain components are very small within the 
limit of elasticity, and the usual first approximation is entirely satis- 
factory. However, if at least one dimension of the body is small com- 
pared to the others, as in the case of plates, shells, and beams, it is 
possible that some of the nine derivatives 0£;/dx; are large, although 
the six components of the strain tensor are small. This takes place in 
some of the problems mentioned before in this section. It appears that 
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the most important task of the general transformation theory of finite 
deformations should be a systematic and exact verification of the 
differential equations for large deflection of beams, plates and shells. 
The present deduction of these equations is based on intuitive as- 
sumptions rather than on an exact application of the general equa- 
tions of elasticity for a three-dimensional continuum. 


5. Stress distribution in plastic and noncoherent masses. Let us 
drop now the assumption that the material obeys Hooke’s law. Two 
types of nonlinear strain-stress relations have practical importance. 
First, in some materials, although there is a unique and reversible 
relation between stress and strain, the linear law does not approximate 
this relation with a sufficient accuracy. Second, we find that beyond 
a certain stress limit the deformation is no longer purely elastic; that 
is, after removing the load, we obtain a so-called permanent set. The 
nonreversible part of the deformation is known as plastic deformation 
or plastic flow; the stress limit beyond which plastic flow is obtained 
is called the yield point. In general if plastic flow takes place the yield 
point is raised, that is, additional plastic flow can only be produced by 
further increase of the load. This effect is called cold hardening of the 
material. Stress distribution in bodies with nonlinear elastic be- 
haviour’ and in materials undergoing cold hardening were treated 
analytically only incidentally. However, the limiting case in which 
cold hardening can be neglected is readily accessible to analytical 
treatment and has been the subject of quite a large number of in- 
genious investigations. 

A material that is capable of plastic deformation without cold hard- 
ening is called perfectly plastic. The experiments show that the yield 
condition of a perfectly plastic material is given by a constant differ- 
ence between the largest and smallest principal stress. This condition 
must be satisfied in the whole domain in which plastic flow takes 
place. Let us restrict ourselves to two-dimensional problems, assum- 
ing that the material flows in planes parallel to the x, y plane. In this 
case the value of the third principal stress, which is parallel to the 
z axis, necessarily lies between the values of the two other principal 
stresses. If we derive the stresses oz, oy, Tz, from the stress function 
F(x, y) by means of the relations 

OF 


=> — 4, = — =7, = — = 


the square of the difference of the principal stresses 0; and a2 is equal 
to 
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(5.2) (o, — o2)? = (r + #)? — 4(rt — 5”). 
Hence, the yield condition can be written in the form 
(5.3) (r + i)? — 4(rt — s*) = K’, 


where K is a characteristic constant for the material. 
Equation (5.3) is a nonlinear partial differential equation of the 


if 


Fic. 7 
From A. Nadai, Plasticity 


second order of the hyperbolic type. The equation of the characteris- 
tics is 


dy\? dy 
(5.4) (r — t) — 4s = 
(5.5) 2% 
dx dx 2s 


Let us denote in an arbitrary point the angle between the tangents to 
the characteristic curves passing through this point and the positive 
x axis by 6; and 62; then according to equation (5.5) 0, and @2 satisfy 
the equation 


r—t 
(5.6) tan 20 = 


2s 


Now the angles ¢; and ¢2 between the directions of the principal 
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stresses o; and o2 and the x axis are the roots of the equation 


oy 2s 
(5:7) tan 26 = = 


Oz — Cy r—t 


From (5.6) and (5.7) we obtain tan 20 tan 26=1. The curves which 
follow the direction of the principal stresses are known as isostatics. 
Hence, the characteristics are a system of orthogonal curves that in- 


From A. Nadai, Plasticity 


tersect the isostatics at 45° and are tangential to the direction of the 
maximum shearing stress. Now the plastic flow consists essentially 
of a gliding along the planes in which the maximum shearing stresses 
occur. Therefore the characteristics are the traces of the gliding sur- 
faces in the x, y plane. Sometimes they are visible on the surface of the 
body subjected to plastic flow without any artifice (Fig. 7); in other 
cases they can be brought out by etching (Fig. 8). 

The geometrical properties of the characteristics given by equation 
(5.5) and the integration of the differential equation (5.3) by means 
of Riemann’s method were studied by H. Hencky, L. Prandtl, 
C. Carathéodory and E. Schmidt. Let us introduce the parameters of 


: 
compression 
dd 
7 
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the characteristics a and B as independent variables, the average nor- 
mal stess p= 3(r+/), and the coefficients A and B of the line element 


(5.8) ds? = Ada* + Bdp? 

as unknown functions. Then we obtain the equations 
3? K 

(5.9) =— 
0a0B 

and 

(5.10) = flo), — —= 


where f(a) and g(8) are arbitrary functions of the parameters a and B. 
To integrate the equations (5.9) and (5.10) we start from two given 
characteristic curves a=0 and B=0 and choose the values of A (a, 0) 
and B(0, 8) along these curves, putting, for example, A (a, 0) = B(0, B) 
=1. Since the characteristic curves are orthogonal, 


where @ is the angle between the x axis and the tangent of a curve 
8=const., the tangent being drawn in the direction of increasing a; 
similarly 


Hence, the functions f(a) and g(8) are known and the system of equa- 
tions (5.10) can be integrated by Riemann’s method. If the network 
of characteristic curves is so determined, equation (5.9) gives the 
value of p. Then the magnitudes of the principal stresses can be cal- 
culated from }(¢:+02) =p and |o,—02| =K. The stresses are directed 
at every point parallel to the bisectors of the characteristic curves. 

The stress distribution in a mass of sand which is in the limiting 
state of equilibrium, that is, just about to collapse, can be determined 
in a similar manner. If the sand is under a large external load and the 
influence of its own weight can be neglected, the stress function 
F(x, y) satisfies the equation 


(5.11) (rt — = 0 


i=f 


where f is the coefficient of friction between the sand particles. This 
equation is identical to the fundamental equation of the theory of 


Be 
1 8B a0 
A da op 
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earth pressure developed by Rankine and can be applied to various 
problems in soil mechanics. However, in many applications gravity 
has to be taken into account. The differential equation of the stress 
function in this case is no longer homogeneous, but contains also 
terms which are linear in 7, ¢ and s. Fortunately, in a few cases of 
practical importance the partial differential equation can be reduced 
to an ordinary equation; such problems have been treated by several 
authors. 

Both in the theory of plastic and noncoherent masses, so far we 
have assumed that the stress distribution is two-dimensional. If we 
try to extend the theory to the three-dimensional case, the difficulty 
arises that the gliding conditions (5.3) and (5.11) involve only two 
principal stresses leaving the third principal stress undetermined. 
There are two ways of attacking the problem. First, one can consider 
the equilibrium of the whole system as a variation problem: the po- 
tential energy of the system has to be a minimum under the restric- 
tion that the yield condition is not violated at any point. Hence, we 
have to solve a variation problem, where the extremals satisfy an in- 
equality. The solution of this variation problem leads to the equation 
(5.3) or (5.11) for the largest and smallest principal stresses and to an 
additiorfal equation for the median principal stress. 

The second method of dealing with the problems of three-dimen- 
sional stress distribution is to introduce certain empirical stress and 
rate of strain relations for the plastic deformation. One can assume, 
for example, proportionality between the stress tensor and the deriva- 
tive of the strain tensor with respect to time. This assumption replaces 
Hooke’s law and furnishes the necessary compatibility equations for 
the determination of the third principal stress. The same two meth- 
ods can be applied to the practically quite important case in which 
one part of the body is in the elastic and another part in the plastic 
state. Both methods furnish a sufficient number of boundary condi- 
tions for the surface separating the elastic and plastic domains. How- 
ever, this surface is not a priori given and, therefore, the solution of 
the problem is in general very difficult. It has been worked out only 
in a few geometrically simple cases. 


6. Hydrodynamics of ideal incompressible fluids. In this and the 
following sections we will consider nonlinear problems occurring in 
fluid mechanics. In §§6 to 8 we assume a nonviscous and incompres- 
sible fluid; in §9 we shall take into account viscosity, and in §10 com- 
pressibility. 

The velocity distribution in an ideal incompressible fluid can be 
determined by the so-called Eulerian equations of motion. Denoting 
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the components of the velocity by u;, the density of the fluid by p, 
the pressure by p and the components of the volume forces per unit 
mass by F;, we have three equations for the dynamic equilibrium: 


i Ou; 1 
(6.1) 


and the kinematic equation of continuity: 


(6.2) 
Ox; 

The equations (6.1) are nonlinear in the u;’s. However, if we assume 
that the flow is without vorticity, that is, the u;’s are derivatives of 
a velocity potential ¢, and that the forces F; are derivatives of a force 
potential V, the three equations (6.1) are satisfied by their first inte- 
gral 

1 dp p 


— — —-+W— — V = const., 


6.3) 
( at 2 OX; Ox; p 


and the equation of continuity becomes 
(6.4) Ad = 0. 


Hence, in the case of irrotational motion the problem of determin- 
ing the velocity distribution is reduced to finding solutions of La- 
place’s equation (6.4) that satisfy certain boundary conditions. These 
boundary conditions are in certain cases linear. For example, if the 
flow takes place between solid boundaries, the component of the 
velocity normal to these boundaries vanishes, and, therefore, 0¢/dn 
=0. Hence, in this case the problem of finding the velocity distribu- 
tion is reduced to a well known linear boundary problem of the po- 
tential theory. 

If the fluid is bounded by a free surface, that is, by another fluid of 
negligible specific gravity which is at rest, the boundary condition has 
to express that the pressure p is constant along the free boundary. 
Hence, we have from equation (6.3) 

1 d¢ 
og 4 dp 


— — — — V = const. 
ot 2 Ox; Ox; 


Assuming that gravity is the only force acting on the fluid, we substi- 
tute V= —gx;, where g is the acceleration of gravity directed parallel 
to the negative x3 axis, and obtain 


t 
a 


1940] NONLINEAR ENGINEERING PROBLEMS 643 


0g 1 
ot 2 Ox; Ox; 


(6.5) + gx3 = const. 


The difficulty of solving flow problems with free boundary arises not 
only from the nonlinear character of the boundary condition (6.5), 
but also from the fact that the boundary itself is unknown. Three- 
dimensional problems of this type have been solved only in a few 
cases; hence, we shall restrict ourselves to two-dimensional flow in the 
%1, x3 plane. We also assume that the flow is stationary, that is, ¢ is 
independent of time and therefore 0¢/0t=0. 

Replacing the coordinates x, and x3 by x and y, we have at a free 
surface the boundary condition 


17 /ad\? 
(6.6) =|($) + &) | + gy = const. 


Let us introduce the complex variable z=x+7y and two analytical 
functions of z, namely, f=¢+72, where @¢ is the velocity potential 
and y the stream function, and w=0++7 log q, where @ is the angle be- 
tween the velocity vector and the x axis, and g the magnitude of the 
velocity-vector. The components of the velocity of flow in the x and y 
direction are equal to u=0¢/dx and v=0¢/dy respectively. Hence, 
u—iv=df/dz and, therefore, w= log (u—iv). Then @=tan—'! v/u and 


09\? 
log gq = log | u — iv| = 4 tog | (*) + (=) | 
Ox oy 


are harmonic functions of x and y and also of @ and yw. Putting 
T= log g, we write the boundary condition (6.6) in the form 


(6.7) Ze" + gy = const. 
or, by differentiating with respect to the arc along the free boundary, 
(6.8) Or oy 0 
— =U. 
Os . Os 


Taking into account that the free boundary is a streamline (y =const.) 
Or Or Or 
q = — 


Os 0d 
On the other hand, dy/ds =sin 6 and, therefore, we obtain from (6.8) 


1 @ 
(6.9) — — (e**) = — gsin#@. 


= 
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We now consider ¢ and y as independent variables; then we have to 
satisfy the boundary (6.9) condition along a horizontal straight line 
w=const. in the ¢, y plane, where A9=0 and Ar=0. We remember 
that @ and 7 are the real and imaginary parts of an analytical function 
of 

The nonlinear boundary condition (6.9) can be linearized in the fol- 
lowing two cases: 

(a) If the effect of gravity is negligible and the fixed boundaries are 
plane surfaces, then we obtain for the free boundary, from (6.9), 
(6.10) 

ag — 

The fixed boundaries being plane surfaces, their traces in the x, y 
plane are straight lines, that is, along such boundaries @ is constant, 
and we have therefore 

(6.11) 

This case was treated in great detail by G. Kirchhoff, H. Helm- 
holtz, Lord Rayleigh and others, and has been applied to many prac- 
tical problems such as the vena contracta, the flow of water over a 
weir, the air resistance of an inclined flat plate, the planing of a flat 
plate over a water surface, and so on. 

(b) [f gravity is taken into account and the motion is only slightly 
different from a uniform parallel flow with the velocity U, we can as- 
sume that @ is small and r= 7 +u/U, where T=log U. Neglecting the 
squares and products of 0, u and v, we have from equation (6.9) 


(6.12) 20 
dg 
or 
(6.13) + =0. 
oy 
Let us assume that Y=0 corresponds to the free surface, Y= — © to 


the infinite depth of the water mass. We also assume that 0, u and v 
disappear at infinite depth. Then 


(6.14) 0 = sin Ap 


satisfies the equation A@=0 and the boundary condition for Y= — «. 
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We still have to satisfy the boundary condition (6.13) for ~=0. Sub- 

stituting (6.14) into equation (6.13), we have 

(6.15) +g=0 

and, therefore, 

(6.16) 0 = sin gb/U?. 

Let us denote the equation of the free surface by 7=7n(x). With the 


free boundary 


flow direction 


wake 


free boundary 


Fic. 9 


approximation used here we can replace ¢ by Ux and dn/dx by 8. 
Hence, we obtain 


xg 
(6.17) = mo COS 
where yo= —0,U?/g. 

If we consider the flow relative to a coordinate system moving from 
right to left with the velocity U, the surface (6.17) represents a wave 
progressing with the velocity U on the surface of a water mass which 
is at rest at infinite depth. The wave length L is equal to 2xU?/g. 

The solution of waves of small amplitude outlined here is the classi- 
cal solution of Airy and Stokes. In the next two sections we drop the 
restrictions that made the linearization of the boundary conditions 
possible and consider the resulting nonlinear problems. 


7. Fluid jets. T. Levi-Civita was the first to extend the theory of 
Kirchhoff and Helmholtz to jets bounded by curved survaces. In a 
paper published in the Rendiconti del Circolo Matematico di Palermo 
in 1907, he treated the problem of the wake behind a blunt cylindrical 
body (Fig. 9). Let us assume that the fluid follows the curved surfaces 


B 
A 
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AB and AC, and that behind the body a dead fluid region, the so- 
called wake, is bounded by free surfaces which extend from B and C 
to infinity. If we consider @ and y as independent variables, and 0 
and r =log g as two harmonic functions of ¢ and y, the boundary con- 
dition for the free surfaces is r=const. To find the boundary condi- 
tion for a fixed curved surface it is convenient to use so-called natural 
equations. For example, we can express the inclination @ as function 
of the length of arc s. Then the curvature d6/ds is also a known func- 
tion of 6, say d0/ds =kK(@), where k is a numerical parameter. Now 
the trace of the fixed surface is a streamline, and, therefore, 


06 00 

—=> — = e— 

Os Os dg 
Hence, we obtain the boundary condition 


06 
(7.1) — = ke*K(8). 


It is seen that for plane surfaces K(@)=0. Then we have the special 
case treated by Kirchhoff, Helmholtz and Lord Rayleigh and equa- 
tion (7.1) is reduced to equation (6.11). 

Levi-Civita, by use of convenient transformations of the variables, 
reduced the problem to the following: 

To find an analytic function 9+ir of a complex variable { =£+7n, 
which is zero for §=0, regular in the half of the unit circle above the 
real axis, and whose real and imaginary parts satisfy the following 
boundary conditions: 

(a) on the circumference of the unit half circle 


do 


f(o)e-*K(8) 


where o is the argument of { on the unit circle and f(a) is a known 
function of 

(b) on the diameter along the real axis r=0. 

The problem can be solved by successive approximations or by ex- 
pansion in a series of powers of the parameter k. 

The theory of Levi-Civita was further developed by S. Brodetzky, 
C. Schmieden, S. Bergmann, V. Valcovici, R. Thirry, L. Rosenhead 
and others. An interesting aspect of this theory is that the point of 
separation, that is, the point at which the free boundary begins, is 
undetermined between certain limiting values. These limiting values 
are determined by the conditions that the free boundaries shall not 
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intersect either the fixed boundaries or each other. The question of 
existence and uniqueness of the solution for a given separation point 
was discussed by A. Weinstein, G. Hamel, H. Weyl and K. Friedrichs. 

If we take gravity into account and the fixed surfaces are curved, 
the boundary conditions are nonlinear both at the free and the fixed 
boundaries. The theory of such “heavy” jets has many important 
practical applications in the design of weirs, spillways, and so on. A 
spillway (Fig. 10) is a steep channel at the downstream face of a dam, 


Fic. 10 


designed for the smooth discharge of large quantities of water. The 
engineer has to design the spillway in such a way that the pressure 
does not sink below the atmospheric pressure at any point on the 
surface where the water is in contact with the spillway. If we com- 
pare the pressure at the points P and Q—that is, on the opposite sides 
of a cross section normal to the average flow direction—we see that 
the height difference between P and Q causes the pressure to be higher 
at Q whereas the centrifugal force causes the pressure to be higher 
at P. Acorrect design requires an approximate balance between these 
two opposite effects. If the pressure at Q is sub-atmospheric, air en- 
ters between the jet and the wall and might cause an unstable jump- 
ing of the jet between two possible flow patterns; one of them is a 
jet with two free boundaries and the other one a jet with one free and 
one fixed boundary. Cases are known in which the instability of the 
jet has caused such violent oscillations of the air that when the spill- 
way was in use windows were rattled within a radius of several miles 
around the damsite. Many engineers have “grappled” with the math- 
ematical solution of the spillway problem; however, the mathematical 
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analysis has not furnished sufficiently dependable information, so 
that spillway design is based mostly on the results of model experi- 
ments. 


8. Waves with finite amplitude. The theory of waves with infini- 
testimally small amplitude sketched in §6 is due to G. B. Airy (1845) 
and G. G. Stokes (1847). Stokes himself, and later Lord Rayleigh, 
T. H. Havelock, A. G. M. Michell and several other authors gave 
further approximations to the exact shape of a wave of finite ampli- 
tude progressing without change (so-called wave of permanent type). 
However, the systematic treatment of the problem, including proofs 
for the existence and convergence of solutions, has been given only 
recently by T. Levi-Civita and his collaborators, D. J. Struik and 
A. Weinstein. 

The analytical problem involved in the theory of waves of finite 
amplitude progressing on the surface of a water mass of infinite depth 
can be formulated in the following way: 

To find all functions 6+i7 of a complex variable {=pe** that are 
holomorphic in the interior of the unit circle |¢| =1, have the value 
zero at ¢=0, and satisfy at the circumference of the unit circle the 
condition 

dr 
(8.1) — = pe-* sin 0 

da 
where p=gL/2rU? is a positive undetermined parameter. The sym- 
bols 6, g, L and U have the meaning given in §6 and r=log g/U. The 
variable ¢ is connected with the complex potential function f=¢+i 
used in §6 by the relation 


(8.2) UL, 
If r and @ are small, we have approximately dt/do = p@ or 
00 
(8.3) p— = 
Op 
The problem is in this case a linear eigenvalue problem; the solution is 
(8.4) 0 = Oop” sin p(o — ao) 


where # is an integer and 6 and oo are arbitrary constants. However, 
it can be shown that physically the solutions corresponding to differ- 
ent integral values of p are all identical. If we choose, for example, 
p=n, this only means that we consider a group of n identical waves 
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with the wave length L as one wave of the wave length nL. Hence, 
it is sufficient to consider the case p=1, which gives the relation 


(8.5) U? = — 


between the velocity of propagation U and the wave length L, in 


troug 
Fie. 11 


accordance with the result given in §6. The factor in equation (8.4) 
and, therefore, the height of the waves in this approximation remain 
undetermined. 

Levi-Civita has shown that if p<1 but 1— 4 is sufficiently small, a 
certain wave shape of definite height/length ratio corresponds to 
every value of p. He also obtained convergent series for the calcula- 
tion of the wave shape. He uses as a first approximation the solution 
of the linear problem for p=1, namely, 


(8.6) 6+ ir = — it = — + yupsing 


where yp is a real constant. Then the solution which satisfies the exact 
boundary condition (8.1) can be written in the form 


(8.7) 6+ i7 = > 


where the w,(¢)’s are polynomials of the mth degree in £. Levi-Civita 
calculated numerically the polynomials w,({) up to »=5. Fig. 11 
shows an example of the wave shape resulting from his computa- 
tions compared with the shape of an infinitesimally small wave of the 
same velocity of propagation. The limiting smallest value of p and 
the corresponding shape of the highest wave has not as yet been ex- 
actly determined. G. Stokes believed that the highest wave has a sharp 
edge at the crest with a vertex angle of 120°. 

D. J. Struik extended the investigations of Levi-Civita to waves of 
finite height in canals of finite depth. A. Weinstein investigated the 
interesting problem of the so-called “solitary wave,” that is, a single 
elevation of the water level travelling with constant velocity and pre- 
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serving the same shape along a canal of finite depth. This phenome- 
non was studied as early as 1844 by J. Scott Russel, who observed 
disturbances of this type generated by the motion of large barges in 


Little Colo. sile gms /litre 
Slope 1% 
Temperature 86 ° 
Rate of flow 750 cc /minute 


YF = 0.0018 
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Courtesy of H. S. Bell, Soil 
Conservation Service 


relatively narrow canals. He followed them several miles on horse- 
back and measured in this way their velocity of propagation. 

Nonlinear problems of fluid motion in rivers and artificial channels 
were treated in the “hydraulic manner” by J. Boussinesq. The “hy- 
draulic manner” of dealing with flow problems consists essentially 
in replacing the partial differential equations for the velocity distribu- 
tion by an ordinary differential equation for the mean velocity aver- 
aged over a cross section. By combining this dynamical equation with 
the equation of continuity, one obtains an ordinary nonlinear differ- 
ential equation for the height of the water level. A great number of 
problems treated in this way can be found, for example, in P. Forch- 
heimer’s book Hydraulik, Teubner, 1930. 

Interesting nonlinear problems related analytically and physically 
to the problem of surface waves arise in the analysis of the flow in 
stratified fluid media. The stability of heavy fluids of different densi- 
ties flowing one above the other with different velocities was investi- 
gated by means of the method of small oscillations by H. Helmholtz; 
V. Bjerknes considered the small perturbations occurring along plane 
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surfaces of discontinuity between fluids of different densities for sev- 
eral combinations of density and velocity distribution. However, the 
exact shape of separating surfaces with large curvature has not yet 
been thoroughly investigated. 


—> 


x/ 3 
YY 
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A case of this type is, for example, the intrusion of a heavy fluid 
mass moving with constant velocity into a mass of lighter fluid, as 
illustrated by the meteorological phenomenon of a progressing cold 
front or by the intrusion of muddy water into a reservoir filled with 
clear water, or salt water into amass of fresh water. Fig. 12 shows the 
shape of the surface separating muddy and clear water. The front of 
muddy water progresses with approximately constant speed and 
keeps its permanent shape. Hence, if we consider the flow relative to 
a coordinate system moving with the heavier fluid, the latter can be 
assumed to be at rest and the lighter fluid to be in stationary motion 
(Fig. 13). The pressure p must be equal on both sides of the discon- 
tinuity surface. Then, if the specific gravity of the heavier fluid is y: 
and that of the lighter fluid y2, the equilibrium condition for the 
heavy fluid yields the condition p/y:+y=const., and Bernoulli’s 
equation applied to the lighter fluid leads to p/y2+y+4?/(2g) =const., 
where y is the height measured from the bottom and q is the velocity 
of the lighter fluid along the discontinuity surface. Taking into ac- 
count that, at the intersection of the discontinuity surface and the 
bottom, y =0 and g=0, we obtain the condition for the discontinuity 
surface: 
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(8.8) — — —— y= 0. 
2g v2 

It is seen that (8.8) is a nonlinear boundary condition of the type of 
equation (6.6). The main interest is directed toward the determina- 
tion of the shape of the discontinuity surface. We can show—using a 
reasoning similar to that employed by G. G. Stokes for the determina- 
tion of the steepest slope occurring in waves of finite height—that if 
frictional effects are neglected the slope of the discontinuity surface at 
the intersection with the bottom must be equal to 7/3. The complex 
velocity potential for the flow of the lighter fluid in the neighborhood 
of the intersection z=x+iy=0 has the form 


(8.9) f = 


where aand C are real, positive constants. Hence, the magnitude of 
the velocity is given by |df/dz| =Ca|z*-'|, and, therefore, q is 
proportional to y*-! and g? to y*‘*-», Since, according to equation 
(8.8), g? is proportional to y, we must have 2a—2=1 or a=§#. Hence, 
f=Csz*!? or z=(f/C)?!*. Since the streamline =0 consists of the bot- 
tom line and the trace of the discontinuity surface, and f=0 corre- 
sponds to the point of intersection, the angle between the two 
branches of the streamline y = 0 is equal to the change in the value of 
the imaginary part of log z= log (f/C), when we pass from ¢ <0 to 
that is, equal to 27/3 =120°. 

In an actual experiment (Fig. 12) the friction on the bottom causes 
a certain distortion of the discontinuity surface; however, if this fric- 
tion is diminished by introducing a thin fluid layer at rest under the 
progressing muddy water, the slope obtained experimentally is in fair 
agreement with the calculated value. It can also be shown analyti- 
cally that the front of the heavy fluid must have a so-called head, 
whose peak is considerably higher than the mean thickness of the 
heavy layer, which is equal to 

U? 
y = — 
Va 2g 

where U is the velocity of progression of the front. However, the 
shape of the discontinuity surface has not yet been determined; it is 
not known whether a mathematical analysis would lead to a surface 
approaching the horizontal level asymptotically or to waves with de- 
creasing or constant amplitudes. 

Also, the nonlinear boundary problems resulting from the analysis 
of the motion of planing surfaces have not yet been attacked seri- 
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ously. The shape of the free water surface behind a planing surface 
has been determined by H. Wagner, but he neglected the effect of 
gravity. The investigators who have taken into account gravity® as- 
sumed small elevations and depressions of the water surface. The gap 
between these two first approximations ought to be filled by further 
investigations. 


9. Viscous fluids. The dynamical equations for the flow of a viscous 
incompressible fluid have the form: 
Ou; Ou; 1 Op 


at * Ox, p ax; 


where 1=1, 2, 3, and v is the kinematic viscosity of the fluid; the con- 
tinuity equation is, in accordance with equation (6.2) 


Ou; 
(9.2) —=0 


Equations (9.1) are known as the equations of Navier and Stokes. 

In the case of viscous fluids we have to assume that at fixed bound- 
aries not only the normal but also the tangential components of the 
velocity vanish. The equations (9.1) and (9.2) are to be integrated 
with these boundary conditions. 

There are a few cases in which this system of nonlinear equations 
tan be exactly solved by known functiorts:” 

(a) In some cases the solutions which satisfy the linear part of the 
equation make the nonlinear terms vanish identically. Thus the prob- 
lem is reduced to the integration of a system of linear equations. For 
example, if only one component of the velocity, for example, the com- 
ponent 1, is different from zero, from equation (9.2) it follows that 
u, is a function of x2 and x; only. Let us consider stationary flow; then 
we have to conclude from equations (9.1) that 0p/dx.=0p/dx3;=0 
and 


(9 3) ) 1 dp 
‘ Ox; p dx, 


The solution of this equation furnishes the velocity distribution in a 
cylindrical pipe. The boundary condition is u;=0 at the wall> 

(b) Certain solutions of the equations for irrotational flow of non- 
viscous fluids lead to exact solutions of the system of equations (9.1) 
and (9.2) if we assume that the streamlines remain unaltered, while 


3 Cf., for example, Lamb’s Hydrodynamics, pp. 375-383. 
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the magnitude of the velocity changes. A typical flow of this type is 
the two-dimensional radial flow between two divergent or convergent 
walls. Writing equations (9.1) and (9.2) in polar coordinates and de- 
noting the radial velocity by u,, we obtain the dynamical equations 
for stationary flow in the form 


du, 1 ap . 
v 


or p on Or? r Or 36? 


(9.4) 


while the continuity equation is solved by 
(9.5) u,r = f(8). 


We introduce f(@) as the unknown function, substituting u,=f(0)/r 
in equations (9.4), and eliminate the pressure p from these two equa- 
tions. Then we obtain for f(@) the ordinary differential equation: 


(9.6) = const., 


which can be integrated by use of elliptic functions. The boundary 
conditions are u,=0 for 6= +a, where a is half of the angle of con- 
vergence or divergence. The interesting feature of this problem is that 
if the flow is convergent, equation (9.6) has only one solution which 
satisfies the boundary conditions at the walls; however, if the flow is 
divergent and the total amount of fluid flowing through a cross sec- 
tion Q=rf*t&u,d0 exceeds a certain value, the solution is no longer 
unique. For example, for the same value of Q, we may obtain a sym- 
metrical and two asymmetrical solutions. The symmetrical solution 
represents a flow in which the fluid shoots with high velocity through 
the center of the diverging channel and back flow takes place at both 
sides of the jet. The two asymmetrical solutions represent flow pat- 
terns with positive velocity near one wall and backflow in the neigh- 
borhood of the other wall. For very large values of Q the sign of the 
velocity changes many times between = —a and =a. This solution 
was found and discussed by the present writer in collaboration with 
one of his Ph.D. candidates in 1913. The young man was killed at the 
beginning of the first world war and the investigation remained un- 
published. Shortly afterward the same solution was found by 
G. Hamel as a special case of a wider class of exact solutions of 
the Navier-Stokes equations. 

Methods for approximate integration of the differential equations 


1 1 Ou, 
r? 00 
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(9.1) are known for two limiting cases: namely, for very large and 
very small values of the parameter rv. 

(A) If v is large, we obtain a first approximation by neglecting the 
quadratic terms, so that we have 


Eliminating p and using equation (9.2), we obtain three linear differ- 
ential equations for 1, u2 and u3. These equations were integrated for 
certain geometrically simple cases by G. G. Stokes. The second ap- 
proximation should be obtained by substituting the first approxima- 
tion in the quadratic terms. It is difficult to find a satisfactory justifi- 
cation for this method of approximation. As a matter of fact, it was 
shown by C. W. Oseen that in the case of a solid moving in a fluid, the 
solution of (9.7) represents a first approximation only in the neighbor- 
hood of the solid, no matter how large the parameter v may be. 

(B) From an engineering point of view the limiting case of small 
viscosity is more important. This topic is known as the theory of the 
laminar boundary layer in the engineering literature. The basic math- 
ematical idea of this theory is the approximate solution of the 
hydrodynamic equations by asymptotic integration. L. Prandtl’s 
fundamental paper On the theory. of fluids with small viscosity was 
presented at the First International Congress for Mathematics in 1904 
in Heidelberg and drew very.little attention from the mathematicians 
attending this congress who were much more concerned with the 
problems of George Cantor’s Mengenlehre and other subjects of pure 
mathematics, on which interest was focused in those days, than in 
the speculations of a young professor of applied mechanics, whose 
presentation perhaps was not very clear from a purely mathematical 
point of view. Several years passed before the method was correctly 
understood. Nowadays it is widely used in the theory of skin friction 
and many other hydrodynamical, aerodynamical and meteorological 
phenomena in which inertia forces and viscous forces are of the same 
order of magnitude. 

Let us consider the following simple problem: An incompressible 
viscous fluid flows along a flat plate that extends along the positive 
x axis from x=0 to x= © (Fig. 14). (We shall use the coordinates 
x, y, 2 instead of x1, x2, x3.) No other boundaries exist; the velocity of 
flow at x= — © is directed parallel to the x axis and is equal to U. 
We assume two-dimensional stationary motion and introduce the 
stream function ¥(x, y) defined by the relations 


+ 
p Ox; 
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which satisfy identically the continuity equation. Then the dynami- 
cal equations yield the following differential equation of fourth order 


for y: 
0. 
(9.8) 


If y=0, the presence of the plate has no influence on the motion and 


J 
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we have u=U or ~o= Uy; the only boundary condition at y=0 is 
v=0. However, if v0, no matter how small v may be, we have the 
additional boundary condition u=0 for y=0. 

To obtain an approximate solution for small values of v, we use the 
method of perturbations, starting from the solution for »=0. We no- 
tice that the perturbation term contains higher derivatives of y than 
the differential equation (9.8) with »=0. It is known that in such a 
case we have to change the independent variable by a substitution 
which contains the parameter v so that we write 


(9.9) = Uy + v(x, 0) 


where n»=yr and a is positive. The exponents a and 8 are to be de- 
termined in such a way that u=0y/dy is of the order of U and the 
largest terms on the left and right sides of equation (9.8) are of the 
same order. The first condition leads to a+8=0. Substituting the 
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expression (9.9) into equation (9.8) and putting 8 = —a, we find that 
the largest terms on the left side are 


Oy Axdy? On 


and 


the largest term on the right side is v 04Y/dy* = v'—**04Y,/dn*. These 
terms are of the same order if a= 4. Hence equation (9.9) becomes 


(9.10) ¥ = + vY,(x, 7) 


where n=y/v'!2. 
With this substitution, equation (9.8) takes the form 


(04 28) th 
On Oxdn? Ax An? On! 


(9.11) 


This equation can be integrated once with respect to 7: 


U + — -— = + f(x). 
On dx On? f(z) 


Since dy,/0n and 0%,/dn? are zero for n= ©, the arbitrary function 
f(x) must be zero. We notice that the substitution ¥(x, 7) =x"/2f(n/x"/?) 
reduces equation (9.12) to an ordinary differential equation for the 
function f where {=7/x"/?. In fact, we obtain 


(9.13) =f". 


This differential equation has been integrated numerically for the 
boundary conditions f=0, f’=—U at and f’=0 at by 
H. Blasius (1908) and more accurately by C. Tépfer (1912). The re- 
gion near the wall where f or the stream function of the perturbation, 
¥i(x, y), is materially different from zero is called the “boundary 
layer.” The order of magnitude of the extension of this region normal 
to the wall is (vx/U)*/?; we say in engineering language that we havea 
boundary layer of the thickness 5, where 4 is of the order of (vx/U)*/?. 
To fix ideas, 5 can be defined by the relation v'/?f°dy,/dndn = U5. 
This value of 6 is called the displacement thickness, and is given by 
6=1.72/(vx/U)?. The velocity distribution in the boundary layer 
is given by u= U+f’. The calculated values of u/U are plotted in 


Ox dy? Ox On? 
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Fig. 15 and compared with experiment. The solution of equation 
(9.13) determines also the skin friction, that is, the shearing stress 
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From Hansen, Aachen Abhandlungen, vol. 8 


at the surface; its value per unit area of the plate is variable with x, 
that is, with the distance from the leading edge of the plate. At a dis- 
tance x, we have 


‘Ou pU? 1 
(9.14) T = pv (=) = 0.66414 ——_——_—_ - 
Oy/ 2 (Ux/v)'!? 


The dimensionless quantity Ux/v is called the Reynolds number of 
the plate extending from x=0 to x=x. 

We have carried through this calculation in the manner of the cal 
culus of perturbations. Prandtl considered the hydrodynamic equa- 
tions in the form given by (9.1) and (9.2) and determined the relative 
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orders of magnitude of the terms in these equations, assuming that 
the velocity component u changes from its value U to zero ir. a small 
range near the wall. Denoting the order of magnitude of this range 
by 6, he concludes that 6 must be of the order v/? and 0u/dy, 02u/dy? 
of the order 1/v'/?, 1/v respectively while 0u/0x is of standard order. 
Then from equations (9.1) it follows that the derivative of the pres- 
sure with respect to a coordinate normal to the wall is of the order 
y'/2, hence negligible in the approximation used in this theory. In this 
way we obtain, for example, in the case of uniform flow along a flat 
plate, the itis equations: 


(0.15) Ou Ou + dv 
Ox dy dy? Ox dy 

The boundary conditions are u=v=0 at the wall and u= U at y= ~, 
where U is the velocity “outside of the boundary layer,” that is, the 
velocity resulting for y=0 from the solution of the equation Ayy>=0 
for nonviscous fluids. 

If we have two-dimensional flow along a wall but the velocity U 
is not constant, equations (9.15) take the form: 

Ou dU Ou dv 


( dx ” ay? Ox dy 


Prandtl has shown that these equations can also be applied to the 
flow along a curved surface, provided the lines x =const. are normal 
to the surface, x is the length. of arc measured at the surface and y 
is the normal distance from the surface. Equations (9.16) furnish a 
fair approximation if K6 and (dK /dx)& are sufficiently small, K the 
curvature of the wall and 6 the thickness of the boundary layer. 

The integration of the system of equations (9.16) has great impor- 
tance for many applications. Various attempts have been made to 
obtain approximative solutions. An excellent review of the methods 
applied is given in Modern Developments in Fluid Dynamics, edited by 
S. Goldstein. Here only a short account will be given: 

(a) S. Goldstein showed that equations (9.16) can be reduced to an 
ordinary differential equation if U(x) has the form const. x* or 
const. e*. The case =0 leads to the case treated above; m=1 that 
is, U=const. x corresponds to a boundary layer starting from a stag- 
nation point of the main flow. This problem has been treated by 
Blasius, Hiemenz, Howarth and Bickley. The equation for arbitrary 
n has been studied by Falkner and Skan and by Hartree. Special 
cases were discussed by R. Mills. 
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(b) Let us assume that U(x) is given in the form of a power series 


(9.17) U = Bix + Box? + Bax? +---. 
Then we put u=0y/dy, v= —dyp/dx and 
(9.18) Fix + Fox? + 


By substituting (9.17) and (9.18) in equations (9.16) and equating 
the terms with the same power of x, we obtain systems of ordinary 
differential equations for the functions F;, F2,--- of y. This method 
was used by Hiemenz, Howarth, Bickley and others. 

(c) We might start from an “initial velocity profile” at x=0, as- 
suming, for example, 

and calculate the “next profile,” that is, u(Ax, y). However, in general, 
that is, for arbitrary values of a1, a2, - - - , we obtain a singularity at 
x=y=O0. If the function U(dU/dx) has the form 


dU 
x 


the function ¥(x, y) is regular only if certain relations between the 
coefficients a1, @,--- and 71,2, are satisfied. Hence, the appli- 
cation of a step-by-step method progressing in the x-direction is 
rather difficult. 

(d) A method of obtaining approximate values of the thickness of 
the boundary layer and the skin friction as function of x has been 
proposed by the present writer and worked out in detail by Pohl- 
hausen. By integrating the first equation of (9.16) between y =0 and 
an arbitrary fixed value of y and using the condition of continuity, 
we obtain 


— u [uv] =v| — —dy. 
( dx 0 4 "I f oy 0 dx 0 2 7 


Taking into account the second equation (9.16), the second term on 
the left side of (9.20) becomes — U(d/dx) fRudx and we obtain 


dx 0 dx 0 dy 0 


Extending the integral to a value y= 6, where U—vw is negligible, that 


= 
E 
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is, over the whole cross section of the boundary layer, we can assume 
that du/dy=0 for y=5 and we have 


(9.23) — u*)dy — — u)dy = (=) . 


This equation, given first by the present writer, is an expression 
for the conservation of the momentum of the fluid involved in the 
boundary layer. Let us assume that the velocity distribution between 
y =0 and y=6 can be expressed by a polynomial of the form 


where a, b and ¢ are functions of x. Then from U—u=0 and du/dy=0 
at y=6 follows: 1—a+b+c=0 and —a+2b+3c=0. For y=0 we 
have, from the first equation of (9.16), v(0?u/dy*)-+ U(dU/dx) =0; 
hence 


dx 


It is seen that a, b and c are expressed by 6 and dU/dx. Therefore, if 
we substitute (9.24) into equation (9.23), we obtain an ordinary non- 
linear differential equation of first order for 6. 

This method was applied with success in many cases and gave fair 
approximations for engineering purposes. It has the disadvantage 
that it is difficult to estimate the degree of approximation. If the ac- 
tual shape of the velocity distribution profile is very different from 
the assumed polynomial, the result can be misleading. 

(e) R. von Mises suggested a transformation of coordinates for 
equation (9.16), introducing x and W as independent variables and 
Z=U?—u? as the unknown function. Then u(0/0x+(v/u)d/dy) has 
to be replaced by u0/dx and 0/dy by ud/dy. We obtain in this way 

(9.25) = mu 

Ox oy? 
where u=(U?—Z?)"/2 and U is a given function of x. If we compare 
this equation with equation (9.12) it is seen that its order is reduced 
from 3 to 2. However, an additional quadrature is necessary to obtain 
the velocity as function of x and y. von Mises suggested a step-by- 
step method to solve equation (9.25). The difficulty caused by a sin- 
gularity at x =y =0 mentioned under item (c) occurs also in this case. 

(f) If we introduce instead of x the variable ¢ = fo Udx (that is, the 


| 
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value of the potential function obtained by solving the equation 
Ay, =0 for the nonviscous case), equation (9.25) can be written in the 
form 


(9.26) =r— 


In the outer portion of the boundary layer, u/U is not very different 
from unity. Then we find that Z satisfies approximately the heat con- 
duction equation: 


The boundary condition is Z = [ U(x) ]? for Y= ©. The solution sat- 
isfying equation (9.27) and this boundary condition might be called 
the “outer solution.” We shall now calculate the “inner solution,” 
that is, approximate values of Z for small values of ¢, and the two 


solutions shall be joined. We choose as the point of junction the value 
of Y=y; for which 


(9.27) 


u 


ayy 
(inflexion point of the velocity profile). Since 0?Z/dy? = —0?(u?) /dy?, 
we obtain from equation (9.26) 
a2(u?) U az 


u =— 
oy? vy 
To obtain a first approximation we will replace 0Z/d¢ by a function 


of @ and u; then we have to solve an ordinary differential equation 
for u, where the coefficients are functions of ¢. For Y=0 


d 


— = 2U — 
dp ’ 


for the point of junction, 0Z/d¢=0. Hence, we replace 0Z/do by 


where u; is the value of u=(—Z?+ U*)"/? obtained from the outer 
solution. 
This method was given by the writer and C. B. Millikan in a joint 


OZ u OZ 
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paper and applied in detail to the case in which U is a linear function 
of x, that is, U= Uo+Uix, where Uy and U; are constants. 

Other more or less laborious methods for obtaining approximate 
solutions of the boundary layer equations have been suggested by 
L. Howarth, Green, Falkner and Skan. The equations for nonstation- 


From Prandtl-Tietjens, Hydro and Aeromechanics, vol. 2 


ary and for three-dimensional motion have been studied by several 
authors and integrated in certain special cases. 

One result of general importance is shown by all these investiga- 
tions: The flow pattern depends largely on the sign of the function 
U(dU/dx). If U(dU/dx) >0, which means that dp/dx <0, the thick- 
ness of the boundary layer decreases or increases slowly, and the solu- 
tion of equations (9.16) can be continued without difficulty for the 
whole range in which U(x) is given. However, if U(dU/dx) <0, that is, 
dp/dx>0O, we soon come to a point at the surface where du/dy=0. 
Beyond this point the streamline Y=0 has two branches; the flow 
separates from the surface and back flow takes place between the 
two branches of Y=0. This point is known as the separation point. 
It can be clearly recognized in the flow pictures shown in Fig. 16 
which represent successive stages of an accelerated flow around a cy- 


LAL 
= 
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lindrical body. Unfortunately, practically all the methods mentioned 
above converge very slowly, if at all, when we approach the separa- 
tion point, and, therefore, it is very difficult to discuss the mathe- 
matical character of the solution of the boundary layer equations in 
the neighborhood of this point. An additional difficulty is caused by 
the fact that by the separation of the flow from the surface vorticity 
is transported into the free stream and, therefore, the assumption of 
irrotational motion outside the boundary is no more valid. Hence, in 
such cases the perturbation method fails. 

Due to this failure of the method we do not get an answer for one 
of the fundamental questions of the hydrodynamics of real fluids, 
that is: What is the flow pattern of a real fluid around a submerged 
body in the limiting case v0? As a matter of fact this problem is 
still not solved. Consider, for example, two-dimensional flow around 
a circular cylinder. We are not able to decide whether the flow pattern 
for v0 approaches the potential flow of a nonviscous fluid or a sta- 
tionary flow pattern consisting of a vortex-free region and a wake 
with continuously distributed vorticity, as suggested by Oseen, ora 
nonstationary flow pattern with concentrated vortex columns of al- 
ternating circulation, a flow pattern treated by the present author. 
It seems that we have here an example in which the analytical meth- 
ods are not sufficient, at least at the present time, to solve a problem 
of purely analytical character.‘ 


10. Compressible fluids. In the previous sections the fluid was as- 
sumed to be incompressible, that is, the equation of continuity or 
conservation of matter was used in the form 

Ou; 


10.1 = 0. 
(10.1) 


If we drop this restriction and denote the variable density of the fluid 
by p, we have 


op 


10.2 
( ot Ox; 


0. 


4 The reader may wonder why the turbulence problem has not been included in 
this review. This omission is due to the fact that the study of turbulent motion has 
not, as yet, led to a clear mathematical formulation of the fundamental laws for the 
mean motion. The problem of transition from laminar to turbulent motion has been 
discussed mathematically by several authors by means of linearized equations with- 
out reaching a satisfactory agreement with experiment. The adequate treatment of 
the nonlinear equations is bound to contribute essentially to the solution of this im- 
portant problem. 
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In this case we have five unknown quantities: the three velocity com- 
ponents, the density, and the pressure. To solve the problem properly 
we have to introduce thermodynamical concepts. We have to add the 
temperature as an additional unknown variable and two equations: 
the equation of state and the equation of the conservation of energy 
(mechanical plus thermal). However, in many cases the problem can 
be simplified by assuming a definite thermodynamic process, for ex- 
ample, adiabatic or isothermal change of state. This amounts to a 
known relation between p and p. Let us assume, for example, an ideal 
fluid undergoing an adiabatic change of state, then we have the rela- 
tion p=const. p*. Introducing the velocity of sound c defined by 
c?=dp/dp, we obtain the following system of equations 
Ou; Ou; c? Op Op 0(pu;) 


(10.3) —+ = - 
ot Ox; p Ox; ot Ox; 


= 0. 


These nonlinear equations can be linearized in the following cases: 
(a) Assume that the velocity components u; and the deviation of p 

from average value po are small. In this case we obtain 
Ou; c? dp Op Ou; 


10.4 =0, 


The velocity of sound can be considered as constant. Eliminating the 
velocity components u;, we obtain for p the wave equation: 


(b) If the velocity components 2 and u3 are small whereas u, might 
be large but differs only slightly from a mean value U, it follows again 
that the deviation of the density p from a mean value po is small. 
Then putting U+uj, we obtain 


Ox, Po Ox; ot ax Ox; 


(10.6) 


The velocity of sound can again be considered constant, so that elimi- 
nating the velocity components u/ we have 
(10.7) A 
ot Ox, 


Let us consider the stationary case, that is, deo/dt=0. Then we have 


1 


666 THEODORE VON KARMAN [August 


U2\ ap 2» 
10.8 
(10.8) ( ) yi 


2 2 2 
Ox; Ox; 


It is evident that the linear differential equation (10.8) is of the ellip- 
tic type if U<c (subsonic flow) and of the hyperbolic type, if U>c 
(supersonic flow). In the first case equation (10.8) is reduced to La- 
place’s equation by the affine transformation 


(10.9) = xe Xe, x3 X3. 


If U>c, equation (10.8) is of the type of the wave equation and can 
be integrated by use of well known methods developed for the solu- 
tion of this equation. 

Equations (10.5) and (10.8) are often used for approximate solu- 
tion of the system of equations (10.3). However, in many cases these 
approximate solutions do not describe the physical picture closely 
enough and we are forced to search for methods of solving the non- 
linear equations. 

Let us consider stationary irrotational two-dimensional flow, using 
the symbols x and y for the coordinates, and u and »v for the velocity 
components. Then denoting the potential function by $(x, y) we ob- 
tain, by putting 0p/dt=0u/dt=0v/dt=0 and eliminating p from 
equations (10.3), 


(20.10) (1 - - 25 
c? ax? c? dy 


In this equation c? is a given function of g?=u?+v?. Namely, accord- 
ing to the definition of the velocity of sound c? is a given function 
of p; on the other hand, it follows from (10.3) that 

2 


p 


Let us apply Legendre’s transformation introducing the velocity 
components u and v as independent variables and X =xu-+yv—@ as 
the unknown function. We obtain 


(10.12) (1-4) +2 +(1-5) = 0 


c?/ dv? c? dudv au? 


If in the whole domain occupied by the fluid u?+v?>c?, equation 
(10.12) is of the hyperbolic type and has fixed real characteristics in 
the u, v plane. Therefore, in this case the integration of equation 
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(10.12) can be achieved by the methods of characteristics. For the 
construction of the characteristics it is convenient to use the polar co- 
ordinates g=(u?+v*)/? and @=tan-'! v/u in the u, v plane. In these 
coordinates the differential equation of the characteristic curves be- 
comes 

1d c 
(*0.12) 

q (q? pb c?) 1/2 

Since c is a function of g only, equation (10.13) is a differential 

equation of the first order for g(@). Its solution consists of two families 


v 
q 
4 
subsonic 
regime 
supersonic” 
regime 
Fic. 17 


of curves: one family corresponds to the positive, the other to the 
negative sign of the square root. If the isothermal process is assumed, 
c=const. and the characteristics are two families of involutes of the 
circle of the radius c covering the whole domain of the u, v plane out- 
side of this circle. If, according to the so-called polytropic law, p is 
proportional to p*, we obtain epicycloids which touch the two circles 
whose radii are equal to the minimum value of c and the maximum 
value of g respectively. The characteristic curves are indicated in 
Fig. 17. If the characteristics of (10.12) in the u, v plane are known, 
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the next step consists of constructing the corresponding cnaracteristic 
curves of (10.10) in the x, y plane. L. Prandtl and A. Busemann have 
worked out graphical methods for this purpose. 

It is a characteristic property of equation (10.10) that in the super- 
sonic range it can have discontinuous solutions. Both density and 
velocity change discontinuously at certain surfaces. Such discontinui- 
ties occurring in supersonic flow are known as shock waves; they also 
can be considered as waves of finite amplitude if we look at the phe- 
nomenon from the standpoint of an observer moving with the fluid. 
B. Riemann was the first who pointed out (in 1860) that in a compres- 
sible fluid waves of finite amplitude can produce discontinuities. 
Later W. J. M. Rankine, H. Hugoniot, J. Hadamard, E. Jouguet, 
G. I. Taylor and others, worked out in detail the theory of the dis- 
continuous solutions. Engineering applications of the theory of shock 
waves were made by A. Stodola, L. Prandtl, and many other authors. 
Applications to the ballistic problem were made by E. Mach, 
C. Cranz, G. I. Taylor, J. W. Maccoll, and P. S. Epstein. 

If in the whole domain u?+v?=gq? <c?, the problem of integration 
of equation (10.10) can be attacked by the method of successive ap- 
proximations. Let us consider the classical example of the two-di- 
mensional flow of a compressible fluid around a circular cylinder. We 
assume that the velocity of the flow at large distance from the ob- 
stacle is parallel to the x axis and has the magnitude U. We denote 
the velocity of sound corresponding to the density and pressure pre- 
vailing in the undisturbed flow by co and the exponent of the poly- 
tropic law by x. Then the velocity of sound at an arbitrary point, 
where the magnitude of the velocity is equal to q, is given by 


2 


Introducing this relation into equation (10.10) and using polar co- 
ordinates r and 6, we obtain 


2 U?\ar ar 80 


where M= U/co is called Mach’s number and 


or r? 


(10.15) 
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We assume that ¢ can be developed in a series of ascending powers 
of M?: 

(10. 16) ¢ = + 

Obviously, we obtain from equation (10.15) 

(10. 17) Ado = 0. 


It is seen that $o(r, 8) is the solution for incompressible fluid. We put 
the radius of the circle equal to unity. Then the boundary condition 
is 0@o/dr=0 for r=1, and we obtain: 


(10.18) oo = (r + 1/r) cos 8. 
We also develop q? in ascending powers of M? by writing 
(10.19) 


Substituting the expressions (10.16) and (10.19) in equation (10.15) 
and equating the terms multiplied by the same power of M, we have 
the following relations for ¢1, de, - - - 


2U? \or or 00 


1 2 
A¢2 = (k — 1)(qgo — 
(10. 20) 


Or Or 06 00 


2 2 
1 (= 1 d¢; 
2U2\ ar or 00 


2 2 
1 1 d¢o 


The boundary conditions are 0¢/dr=0 for r=1 and 0¢,/0r =0¢2/0r 
=0 for r=. 

The method sketched here is essentially identical to that used by 
Lord Rayleigh, O. Janzen, L. Poggi, C. Kaplan, I. Ismai, K. Tamada 
and Y. Saito. The calculations of the Japanese authors are the most 
reliable ones; they corrected some numerical errors which occurred in 
previous publications. G. I. Taylor worked out an experimental 
method for solving the same problem by successive approximations. 
He employed an electric device analogous to that used for the experi- 
mental solution of Laplace’s equation. 

Another method for solving problems of subsonic flow in compress- 


ame 
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ible fluids was suggested by P. Molenbroek and A. Chaplygin. These 
authors use a transformation of independent variables which is analo- 
gous to the Legendre transformation mentioned before in this section. 
However, instead of using the potential $(x, y) and its Legendre 
transform X(u, v) as unknown functions they use $(x, y) and the 
stream function ¥(x, y) which satisfies identically the continuity 
equation for stationary flow: 


O(pu)  A(pv) 

Ox oy 
We write pu=d/dy and pu= —dy/dx. Then we have 
(10.22) dp = udx + ody, dy = — prvdx + pudy. 


(10.21) 


We resolve these equations for dx and dy and introduce as independ- 
ent variables polar coordinates in the u, v plane, namely, g = (u?+-0?)'/2 
and @=tan~'! v/u. Substituting 


oy oy oy 
d@é = —dqg+—ds, d¥=—dq+—d, »v=qsin# 


into equations (10.22) and resolving them for dx and dy, we obtain 
cos sin 6 cos sin 0 
ax =( + ( a, 
q 99 9% q 96 9 


sin @ cos@ dy sin df cos@ dy 
q 99 q 90 


(1.23) 


The relations 
a (2) a (%) a (2) a (2) 
a0\aq/ aq\ae/’ a0\aq/ aq 
furnish the following two relations between the derivatives 0¢/0q, 


06/06 and df /dq, dy /00: 


. 
— sin 6— + = — — —— 


sin 0 
oy cos@ sin (+2 2)% 


a a dq) 30’ 
(10.24) q qP q 
sin 0 @ cos 6 
sp 0q q 90 qp p | 


Multiplying the first equation of (10.24) by cos @ and wie second 
by sin 0, we have by subtraction 


| 
| 
| 
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(10.25) 


Multiplying the first equation by sin @ and the second by cos @, we ob- 
tain by addition 


1 d 
(10. 26) —(1+4 
4p p dq/ 06 


From equation (10.11) follows (q/p)dp/dq= — q?/c?; hence, equation 
(10.26) becomes 


0 1 Ne 
(10.27) 
qp 
Eliminating, for example, ¢ from equations (10.25) and (10.27), we 
obtain the following linear differential equation for y: 


p ey 
1 — g?/c? dg\p aq 06? 


where c is given as function of g by means of the relations c?=dp/dp, 
dp/p+qdq=0 and p=f(p). The last relation depends on the law of 
the change of state of the fluid. 

Chaplygin has shown that the differential equation (10.29) can be 
reduced to Laplace’s equation by an appropriate choice of the rela- 
tion p=f(p). Laplace’s equation for the polar coordinates g and @ has 
the form 


oy ~) 
( ) (« 
Equation (10.29) is reduced to the form (10.29) if 


d 1 pd 1 — g*/c? 
(10.30) 
q aq qp 


dq 


where a is an arbitrary constant. Hence, we have to choose p=f(p) 
such that 
p> = a°(1 — g?/c*) 


or 


(10.31) g=c— 


p og q 
= 
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Since d(q?)/dp= —2c?/p, we obtain by differentiation of (10.31) the 
relation 


(10.32) 
From equation (10.32) follows c?=dp/dp = B/p? and 


B 
(10.33) pul 

p 
Introducing the specific volume V=1/p, we can write equation 
(10.33) in the form 


(10.34) Pp — po = const. (Vo — V). 


It is seen that we obtain Laplace’s equation for ¥ (and also for ¢) if 
the polytropic curve in the p, V plane is replaced by a straight line. 

It is obvious that if the density change is small the adiabatic or 
other polytropic curves can be approximated by straight lines. 
A. Chaplygin, D. Riabouchinsky, B. Demtschenko, J. Caius and 
A. Busemann applied this method of approximation to the slow mo- 
tion of a compressible fluid. It seems to the writer that an application 
to the flow of a compressible fluid with a velocity which is comparable 
to the velocity of sound is of greater practical interest. Let us assume 
that the undisturbed velocity of the fluid is U (where U<c) and the 
density corresponding to the velocity is pp=1/Vo. If we assume that 
the change of the density is sufficiently small so that the pressure can 
be expanded as function of V— Vo» and approximated by a linear. 
term, it is possible to determine possible flow patterns, velocity and 
pressure distributions in compressible fluids by use of the methods 
of conformal transformation. By means of successive approximations, 
which require considerably less labor than Rayleigh’s method men- 
tioned before in this section, it is also possible to determine the veloc- 
ity and pressure distribution around a given body. H. S. Tsien worked 
out this application of Chaplygin’s method upon suggestion of the 
present writer. 

The problem of flow of a compressible fluid around a cylindrical 
body or a body of axial symmetry is of great practical importance in 
view of the increasing speed of aircraft, which already have flown 
faster than 2/3 of the velocity of sound. Experimental evidence seems 
to show that if the velocity at any point of the wing surface reaches 
the velocity of sound, a shock wave appears. However, G. I. Taylor 
calculated numerically with great patience an analogous flow prob- 
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lem, namely, the flow of a compressible fluid through a nozzle, whose 
cross section decreases to a minimum value and again increases. He 
obtained a flow pattern, showing near the walls a velocity greater 
than that of sound, but no discontinuity. Another important case is 


Fic. 18a. Photograph of a bullet in flight. From G. I. Taylor and 
J. W. Maccoll, Aerodynamic Theory, vol. 3 


the motion of a blunt-shaped body at supersonic velocity. Observa- 
tions show that a curved discontinuity surface is created in front of 
the body; between this surface and the blunt nose of the body the 
flow is subsonic and the velocity becomes zero at the nose of the body. 
The same phenomenon occurs also in the case of a body with a 
pointed nose or leading edge at velocities slightly higher than the 
velocity of sound (Fig. 18a), whereas if the speed is sufficiently large, 
the shock wave starts from the pointed edge, or in the case of a body 
of rotation from the apex (Fig. 18b). This problem is important in 
ballistics; for example, the so-called terminal velocity of aerial bombs 
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dropped from large heights is in this range. Unfortunately, no method 
as yet is available for solving such mixed subsonic-supersonic prob- 


Fic. 18b. Photograph of a conical headed projectile in flight. 
Photograph by T. Harris, British Ordnance Committee 


lems. The difficulty is that the differential equation for the potential 
or the stream function is of the elliptic type in one part of the domain 
and of the hyperbolic type in the other. It seems to the writer that 


wit 
Sin 
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the invention of methods for the solution of such mixed problems 
would be an achievement both from a practical and a mathematical 
point of view. 

An attempt was made in this paper to show the applications of 
analytical methods available for the solution of certain nonlinear 
problems in which the engineer is interested. Some gaps were shown 
and frontiers indicated beyond which the safe guidance of the mathe- 
matical analysis is for the time being lacking. I would consider this 
lecture a great success if it would induce mathematicians to volunteer 
for pioneering work in pushing those frontiers “farther west.” 
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SULLE TRASFORMAZIONI CREMONIANE CHE 
POSSEGGONO PER CURVA DI PUNTI UNITI 
UNA SESTICA CON DIECI PUNTI DOPPI 


GIUSEPPE POMPILJ 


In una precedente memoria Sulla trasformazioni Cremoniane del 
piano che posseggono una curva di punti uniti' ho messo in evidenza 
che tra i tipi possibili di siffatte trasformazioni si presentano quelle 
che hanno per curva di punti uniti una sestica con 10 punti doppi; 
per provare l’effettiva esistenza di tali trasformazioni ne ho dato un 
esempio (num. 21) che perd é risultato essere insussistente perché, 
come ha recentemente osservato A. B. Coble,? la trasformazione che 
si ottiene non é altro che l’identita. Percid do qui un altro esempio. 

La generica sestica razionale, con dieci punti doppi Pi, Pe, - - - , P:, 
A, B, C di cui gli ultimi tre nei vertici del triangolo fondamentale 
A(0, 1, 0), B(1, 0, 0), C(O, 0, 1), ha equazione parametrica: 


t(t — a1)(t — a2)(t — yi)(t — v2), 

t(t — Bi)(t — Be)(t — — v4), 

(t — a1)(t — ae)(t — Bi)(t — Be)(t — ys)(t — ve). 
Imponiamo a questa sestica di soddisfare alla seguente condizione: 

le due involuzioni individuate dalle coppie che cadono nei punti doppi 

A, Ce B, C devono essere permutabili. Questo porta che sia a1a2 = — BiB. 

o meglio B2= —a,a2/B;; cioé si ha un solo legame tra i coefficienti e 

la condizione, come é naturale, risulta perfettamente simmetrica ris- 

petto ai punti doppi A e B. In questo caso I|’equazione della nostra 

sestica S diventa: 

— a)(t — ae)(t — yi)(t — 

t(t — Bi)(¢t + arae/Bi)(t — ys)(t — va), 

(¢ — — a2)(t — Bi)(t + — ys)(t — v6), 


x 


y 


x 


y 


ed un semplice esame di questa equazione parametrica ci assicura che 
nell’intorno dei tre punti doppi A, B, C non esistono altri punti 
multipli di S. 

Ma ancora con una sola condizione si pud imporre che la sestica 
abbia: 


1 Rendiconti del Seminario Matematico della Universita di Roma, (4), vol. 2 
(1938). 
2 This Bulletin, vol. 45 (1939). 
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I. un punto triplo in luogo di tre punti doppi, o 
II. un tacnodo, o 
III. tre punti doppi su una retta, o 
IV. sei punti doppi su una conica, 
gli altri casi, come vedremo, non ci interessano perché oltre che 
P,, ---, Pz fanno intervenire il punto doppio C ovvero, per la sim- 
metria, A e B contemporaneamente. 

Sorge cosi il dubbio che la condizione imposta equivalga a qual- 
cuna di quelle sopra elencate. 

Il primo caso si esclude subito. Supponiamo infatti che la nostra 
sestica abbia un punto triplo P in corrispondenza ai valori 1, 2, 73 
del parametro ¢; in tal caso le rette per P segano ulteriomente S 
secondo una g; che sul parametro ¢ é espressa da: 


(t¢ — a)(t — ae)(t — 51) + ACE — Bi)(t + — 52) =0 ; 


siccome poi questa g; deve contenere anche un gruppo del tipo 
(0, ©, 63) segue che 56,= — 6 (=85). Si osservi ora che le rette per 
A e B segano ulteriormente S secondo due gj espresse, sul parametro 
t dalle equazioni: 


— yi)(t — v2) + A(t — Bi)(t + — ys)(t — = 0 
t(t — ys)(t — ys) + w(t — an)(t — ae)(t — — 0) = 0 


e queste due g; devono possedere rispettivamente due gruppi del tipo 
(5, m, 72, 73), (— 45, m1, 2, 93) il che perd porta un’ulteriore condizione 
sui coefficienti; possiamo quindi concludere che in generale la nostra 
sestica possiede 10 punti doppi e hon un punto triplo e sette punti 
doppi. 

Gli altri tre casi non sono riuscito ad escluderli. Preciserd percid 
che cosa diventano nei vari casi le trasformazioni involutorie del 
quarto tipo di Bertini* che hanno per fondamentali gli otto punti 
Pz, B). 

Se non si presenta nessuno dei casi del numero precedente, ovvero 
se la sestica S ha un tacnodo si hanno due trasformazioni del 17° or- 


dine con i punti fondamentali: (P%, - - - , P§, A®), (P%, ---, P%, B®). 
Se P;, P2, Ps sono. in linea retta si hanno due trasformazioni del 16° 
ordine con i punti fondamentali: (P?, P3, P3, P%,---, P%, 
Se P,,---, Ps sono su una conica si hanno due trasformazioni 
del 13° ordine con i punti fondamentali: (P{,---, P§, P%, A®), 


3V. Snyder, The involutorial birational transformation of the plane, of order 17, 
American Journal of Mathematics, vol. 33 (1911). 
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Cioé indicando con b, la trasformazione del quarto tipo di Bertini 
che ha i punti fondamentali (Pi, --- , P:, A) e con bg quella con i 
punti fondamentali (P;, - - - , Pz, B) possiamo dire che le due tras- 
formazioni b, e bg hanno rispettivamente i punti A e B come fonda- 
mentali e mutano i loro intorni in sestiche che passano doppiamente 
per P;, -- - ,.P;mentre hanno un punto triplol’una in A el’altra in B. 

Data la natura della sestica S le due trasformazioni involutorie 
bs e bg sono permutabili su essa mentre, come dimostreremo, non 
lo sono sul piano; se infatti fosse )4bg =bgb, siccome B é unito per b, 
mentre é mutato da bz in una sestica irriducibile Dg (P?, - - - , P?, B*) 
seguirebbe che questa sestica D, sarebbe unita per b, in contraddizione 
con il lemma VI (num. 10) della mia memoria precedentemente 
citata. Da cid segue che la trasformazione Q =b4bgbabs non é identica 
ed ha la sestica S come luogo di punti uniti. 

Questa trasformazione Q non ha perd tutti e dieci i punti doppi di S 
come fondamentali; per averne una che soddisfi anche a questo re- 
quisito, basta prendere una trasformazione cremoniana T la quale 
oltre che mutare in sé la sestica abbia per fondamentali quei punti 
doppi di S che non lo sono per Q, e quindi costruire la: R= T—*QTQ-". 

R. UNIVERSITA, 

Rome, ITALY 


THE SOLUTION OF BOUNDARY VALUE PROBLEMS BY 
A DOUBLE LAPLACE TRANSFORMATION 


J. C. JAEGER 


1. Introduction. The Laplace transformation has most frequently 
been used to transform a linear partial differential equation with two 
independent variables and constant coefficients into an ordinary dif- 
ferential equation, the “subsidiary equation,” from the solution of 
which that of the original equation is deduced. If it is applied to an 
equation with more than two independent variables, the subsidiary 
equation is also a partial differential equation and is usually solved 
by classical methods. The object of this note is to point out that 
partial differential equations in which the range of two (or more) 
independent variables is (0, ©) may easily be handled by simultane- 
ous Laplace transformations in these variables. The point of view is 
that of previous papers' in which the Laplace transformation method 
has been regarded as purely formal, and the solution as subject to 
verification. The method is related to that of Doetsch? for equations 
of elliptic type but assumes no theoretical basis. In §3 a simple two 
variable problem is discussed to illustrate the method and in §§4 and 
5 are given two new three variable problems of a type to which it is 
well adapted. 


2. The general method. We consider a function v(x, ¢) in the range 
x>0, t>0. 
Let 


Ov 


a(t) = = [= 


Ox. z=0 


The Laplace transform with respect to ¢ will be denoted by a “bar,” 
thus 


Jie) 


and that with respect to x by a capital letter, thus 


1Carslaw and Jaeger, Proceedings of the Cambridge Philosophical Society, 
vol. 35 (1939), p. 394; Proceedings of the London Mathematical Society, in press; 
this Bulletin, vol. 45 (1939), p. 407. 

2 Doetsch, Theorie und Anwendungen der Laplace-Transformation, Berlin, 1937, 
chap. 22. 
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G(p’) -f e~?’*9(x)dx, 
0 
and the double Laplace transform by both: 
Vip’, p) = f f PtP’ 2y( x, t)dxdt; 
0 0 


p and p’ are supposed to have real parts sufficiently large to ensure 
convergence. This notation indicates compactly the variables oc- 
curring; if there is another independent variable y, it will occur 
throughout. 

Then subject to suitable restrictions on v we have 


(1) fie dxdt = pV — Up, 
(2) dxdt = 9°V — 9U. — 
(3) f f dxdt = p'V — wp, 
(4) f J dxdt = — — wh. 


From a given differential equation and boundary conditions we ob- 
tain a subsidiary equation and boundary conditions by multiplying 
by exp [—pt—p’x], integrating with respect to x and ¢ from 0 to , 
and using (1) to (4). Solving this equation we obtain V(p’, p). 

To derive v(x, t) from V(p’, ») we assume an inversion theorem 
which may be derived formally from Fourier’s integral theorem in 
several variables: 

1 yt+ix ico 
(6) af na, 
4x? y—ico 
where >arg A> —z, and >arg > —7, provided that: 
(6) V(X’, d) is bounded in some half-planes R(\’) > a’, R(A) > a, 
and y >a, y' >a’. 


Other conditions in addition to (6) will of course be necessary for 
the truth of the inversion theorem (5), but neither these nor the as- 
sumptions involved in (1) to (4) need be discussed since the whole 
process of determining v(x, t) as a double contour integral of type (5) 
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is regarded as purely formal. It has then to be verified that this solu- 
tion does satisfy the differential equation and boundary conditions. 
This is done in the case of the equation of conduction of heat by 
transforming the paths of integration in (5) to paths beginning in the 
third and ending in the second quadrant, and the verification is then 
performed on these new integrals. The whole process follows the lines 
of the one-variable case* and will not be given here. The contour 
integrals of type (5) are evaluated by a deformation of the contour 
precisely as in the one-variable case, this will be done without com- 
ment, the method and justification having been given in the papers 
referred to. 


3. Heat conduction in the region x>0; the initial temperature 
zero; the end x =0 kept at unit temperature for ‘>0. We have to solve 


with4 

(9) w(t)=1, w,i(#) unknown. 


By (1) and (4) the subsidiary equation is 


x p 


and so 
+ pir(p) 
V(p", p) 
For V(\’, A) to satisfy (6) we must have® #,(A) = — (Ax)~/? and thus 
= 1 
x) = 


+ 


Hence by the inversion theorem (5) 


3 Carslaw and Jaeger, loc. cit. 
4 Boundary conditions stated shortly as in (8) and (9) are to be understood as: 


(8’) lim v(x, t) = 0, for fixed x > 0, 

(9’) lim v(x, t) = 1, for fixed t > 0. 
0 


x 
5 From this result the flux of heat at the origin may be obtained directly. 


v(x, 


1 ico 


2 


gil2 


(10) 


070 
Ox? ay? 

with 
x) = 0, 


t) = 0, 


t) = 0, 
vy, x, 4) = 1, 
The subsidiary equation is 
aV 
dy’ 


(11) 


where 


j. 
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yt ico 


4. Heat conduction in the region x >0, 0<y<b; y=b maintained 
at unit temperature for ‘>0; y=0 and x=0 maintained at zero for 
t>0; the initial temperature zero. We have to solve 


= 0, t>0,%>0,0<y<8, 


x>0,0<y<b, 
t>0,0<y<b, 


w:(y, 4) unknown, 


when y = 0, x > 0,#>0, 
when y= 56, %>0,¢>0. 


+ q°V w(y), 


= p'*—p/x. This is to be solved with V =1/pp’ when 
and V=0 when y=0. The solution is 


Now sin b(A\’*—X/x)¥?, qua 
is, 
(—)"r 


+ 


(13) 


+ n=0, 1,--- 
V(y, \) must vanish at \’ = 


b nry’ 
+f 
0 


sin gy sin gy 
V(y, = f sin g(b — 
(12) pp’ singb qsingbJ, 
sin 
q sin qb 0 


function of \’, has zeros at \’= 
, so to satisfy (6) the numerator of 
+ n=1, 2,---, that 


tir(y’)dy’ = 0, 
n= 1,2,--- 
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The residue of e’*V(y, d), qua function of at the pole 
= is 
exp [— x(m?x?/b? + (—)*"r 
sin —— 
-f sin 
0 b 


Using (13) in this gives for the result of the \’-integration in (5) 
sinh 2xnn(—)” 
d sinh b(A/x)!/? A(kn?x? + db?) 


- sin € 


where the first term comes from the pole \’ =0. 
For the \-integration we require 


1 prtio exp — + 


A(A/Kk + a?) 
1 1 ee du 
= — — — =f sin yx 
u(u? + a?) 
and 


Using these results we have finally 


n 
vy, x, t) = gore sin 


n=l nN 
_ 
(14) + — n(—)* sin 
b? n=1 
du. 
0 u(u? + n?x?/b?) 


Clearly problems in which the temperature of y=) is a function of 
x or t may be dealt with in the same way provided the function satis- 
fies fairly wide conditions. 


5. Heat conduction in the region x>0, 0<y<b; the initial tem- 


ry 
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perature unity; the surface maintained at zero for :>0. We have to 
solve 


with 
u(y, x) = 1, <9 <6, 
wo(y, 4) = 0, t>0,0<y<b, 
unknown, 
v(y, x, 2) =O when y=O and y=5d for x>0 and ¢>0. 
The subsidiary equation is 
av 1 


y— ’ 
dy? q i(y) 


where g? = p’?— p/k. 
This is to be solved with V =0, when y=0 and y=b. The solution is 


sin gb — sin gy — sin g(b — y) 


(y, P) sin ob 


sin gy Ric? 
q sin gbJ, 
sin g(b — y 
2 sin qy'w(y’)dy’. 
q sin gb 0 


The zeros of sin b(\’?—X/x)"/?, qua function of \’, are at 
= +(n*x?/b+X/x)/?, n=0, 1, 2, - - - , and thus, in order that (6) may 
be satisfied, we must have 


1, 2, 


The residue of e’*V(y, \’, X) at the pole \’ = — (n?x?/b?+2/k)"/? is 


b[1 (—)*] ny’ \ 
+f sin b tly \dy 


Using (16) in this the result of the \’-integration in (5) is 


| 
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4b? = 1 
(2r + 1)[x(2r + 1)2x? + 
[ + 1)?x? (2r + 1)ry 
-exzp | — + — 
b? K b 
sinh b(\/x)'/? — sinh y(A/x)!/? — sinh (b — y)(A/x)!/? 
sinh b(A/x)"/2 


where the last term has been derived from the pole \’ =0; there are 
no poles at A’ = + (A/«)"/?. 
Carrying out the \-integration in the usual way, we have finally 
sin (2r + 1)ry/b 
(2r + 1) 


8 
wy, x, t) 3; 
(17) 


© sin ux 
f exp [— xt(u? + (2r + 1)?x?/b?) 
0 


The same method applies to cases in which the initial temperature 
uo(y, x) is a function of x and y under fairly general conditions on the 
function. 

As an example let uo(y, x) =x, for x>0, 0<y<b. 

The only change in V(y, p’, p) is that the first term of (15) is re- 
placed by 


— sin bg + sin gy + sin g(b — y) 
xp’?q? sin bq 


Proceeding in the same way, we find that the only pole of V(y, \’, A), 
qua function of ’, is \’ =0, where the residue is 


sinh b(A/x)!/? — sinh y(A/x)!/? — sinh (6 — y)(A/x)!/2 
X sinh 


Performing the \-integration of (5), we obtain finally 


42 1 (2r + 1)ry 
18 oy, x, = — em 
(18) oy, #, = — 


THE UNIVERSITY OF TASMANIA 


= 


THE MINIMAL NUMBERS OF BINARY FORMS' 
RUFUS OLDENBURGER AND ARTHUR PORGES 


1. Introduction. One of us proved that for certain fields K a 
form F of degree m can be written as a linear combination of mth 
powers of linear forms. Such a combination is termed a representation 
of F and the least possible number of terms in any such representation 
is called the minimal number of F with respect to K. The minimal 
number depends on both F and K. For fields K with characteristic 
greater than n, and binary forms F of degree n, it has been proved? 


that the minimal number ranges over at least 1, 2,---, m, and at 
most 1, 2, - - - , #+1, but the exact range was not determined. In the 
present paper the authors prove that the range is precisely 1,2, - - -,m. 


2. Preliminary lemmas. In what follows we use identity of poly- 
nomials in the usual sense, namely polynomials P and Q are identical 
if the coefficients of P equal the corresponding coefficients of Q. 

Since the order of a field K is greater than m if the characteristic 
of K is greater than m, we have the following lemma. 


LemMaA 1. For a field K with characteristic greater than m a polyno- 
mial P of degree m is equal to a polynomial Q for all values of the vari- 
ables if and only if P and Q are identical. 


An immediate consequence of Lemma 1 is the following lemma. 


LemMA 2. For a field K with characteristic greater than m, a polyno- 
mial P of degree m not identically zero is different from zero for at least 
one set of values of the variables. 


Lema 3. Let K be a field with characteristic greater than m. Let A be 
the determinant 


1 
(1) am be 
G1 


of order m+1, m=1, with elements in K, and suppose that the b’s are 
not all zero. The determinant A is not identically zero in the a’s. 


1 Presented to the Society, April 13, 1940. 
2 R. Oldenburger, Polynomials in several variables, Annals of Mathematics, (2), 
vol. 41 (1940), no. 3, pp. 694-710. 
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Lemma 3 is valid when m=1. For m2=2 we have the following se- 
quence of equalities: 


1 0 tee 0 by 
(d2— 41) ---(dm— a) be 
a, (a2 
---0 bi 
(be — bya) 
= Ge +--+ Om (bs — bear) 
Thus 
(2) A = (a2 — a1) (@m — a1)(M — aN), 
where 
1 be 1 | bi 
2 "0s a2 
M = as On bg. 1; N= as On bs 


The lemma now follows by induction. 

It is to be remarked that if the characteristic of a field K is greater 
than m, each binary form of degree m with coefficients in K can be 
written as 


(3) bix™ + mbex™ ly +--+ + 
where the coefficient of b;x™-‘+1y‘-! is the binomial coefficient of 
x™—itlyi—! in the expansion of (x+)”. 

3. The range of the minimal number. We proceed with the follow- 
ing theorem in which J(m) denotes the integer 
(4) 3m|m(m — 1) + 2]. 

THEOREM 1. For a field K with characteristic greater than m, and 


order at least I(m), the minimal number of a binary form F does not 
exceed m. 
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We write F as in (4). We exclude the trivial case F=0. If 5:+0 
while b;=0 for i>1, the form F is simply b,x”, whence the minimal 
number of F is 1; similarly, if bn4:#0 while 6;=0 for i<m-+1. If 
m=2 while }; and bas, are not zero, and if further 6;=0 for 
i=2, 3,---, m, the form F is nonsingular in the sense of Olden- 
burger,* and has minimal number 2. 

If F is of the first degree, the minimal number of F is clearly 1, 
while if F is quadratic, the minimal number is identical‘ with the 
rank of F and does not exceed 2. 

It remains to consider forms F of degree at least 3 for which 
be, - - - , bm are not all zero. It is no restriction to assume that bm410. 
For if bn» ,:=0, the form F can be transformed nonsingularly into a 
form with this property, as is clear from the following argument. We 
write F= F(x, y). Since F£0, there exist values a and b of x and y 
respectively such that F(a, b) #0. We make the transformation 


x=x2'+ ay’, y = + by’, 


on F to obtain a form F’(x’, y’). Evidently, F’(0, 1) = F(a, b). 
We consider the following equality: 


(5) + mbex™1y + + = D(x + 


i=1 
The equality (5) is identically satisfied if the \’s and a’s are chosen 
so that the following system of linear equations is valid: 
At Ast---+ 
2 2 2 
(6) + + --- + = 
We shall prove that the a’s can be chosen so that the rank of the 
matrix of coefficients of the \’s in (6) is m, while the rank of the aug- 
mented matrix is also m. The determinant of the augmented matrix 
of (6) is A. We write A as in (2). 


In the expansion of M, the terms of lowest degree a in the a’s are 
those which have bn; as coefficient. The degree a is explicitly 


3 R. Oldenburger, Rational equivalence of a form to a sum of pth powers, Transac- 
tions of this Society, vol. 44 (1938), pp. 219-249. 
4M. Bécher, Introduction to Higher Algebra, p. 135. 
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4(m—1)(m—2). The terms of lowest degree, in the expansion of the 
polynomial a;N for any 7, are at least of degree a+1. Thus the terms 
of degree a in the polynomial M—a;N [i in the range 2, 3,---, m| 
arise from M only, whence M—a;N is not identically zero in the a’s. 
Thus the polynomial 


(7) (M — a2.N)--- (M — anN)N 


of degree at most J(m) —1 is not identically zero in the a’s. Evidently, 
N+=0 implies that a2,---, Gm are distinct. Choose a set of values 
of d2,---, @ such that the polynomial (7) is different from zero. 
For this choice of the a’s the polynomial M —a,N is linear in a, and 
M-—a,N=0 has a solution for a, distinct from a2,---, @m. Thus 
there exist mutually distinct values of a;, de, - - - , dm, such that A=0. 
Since the a’s are mutually distinct, the matrix of coefficients of the 
\’s in (6) has rank m. It follows that the augmented matrix has the 
same rank. It is well known that a system of linear equations hasa 
solution if and only if the matrix of coefficients and the augmented 
matrix have the same rank.5 Theorem 1 is now proved. 

The following theorem was proved by Oldenburger.® 


THEOREM 2. For a field with characteristic greater than m, the mini- 
mal number of x™—‘y is m. 


If the sum \,L7 (i =1, 2, - - - , m) is chosen to be a minimal repre- 
sentation of x"~1y, the sum 


MLt p sm, 


is minimal. Thus for each p in the range 1, 2, - - - , m there is a binary 
form of degree m with minimal number p. Applying Theorem 1 we 
have arrived at the following result. 


THEOREM 3. Let a field K be given as in Theorem 1. For the field K 
the range of the minimal numbers of binary forms of degree m is 
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5M. Bécher, Introduction to Higher Algebra, p. 46. 
6 See first reference to Oldenburger above. 
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A REMARK CONCERNING SUFFICIENCY THEOREMS 
FOR THE PROBLEM OF BOLZA' 


E. J. MCSHANE 


In recent years there have appeared several sufficiency theorems 
for minima in Bolza problems in which there is no assumption that 
the curve in question is normal. However, it is assumed in all these 
theorems that the multiplier rule, the transversality conditions and 
the strengthened conditions of Weierstrass, Clebsch and Jacobi hold 
for a set of multipliers with Ay) >0; whereas, among the necessary con- 
ditions the multiplier rule and the Weierstrass and Clebsch conditions 
have been shown? to hold for multipliers with A» 20. 

In the present note we close this slight gap by showing that the 
sufficiency theorems hold if we assume Xo non-negative instead of 
positive. The only case needing discussion is \y»=0. In this case we 
show that the curve satisfying the hypothesis is isolated; no neigh- 
boring curve satisfies the side-equations and end-conditions. The 
curve is thus a minimizing curve in a trivial sense. Less trivially, we 
show’ that the curve is a normal proper minimizing curve for a Mayer 
problem in which the end-conditions consist of a subset of the original 
end-conditions and the function to be minimized is a linear combina- 
tion of the original end-functions. 

The problem is that of minimizing a functional 


(1) = + f "fle, 12), 


in the class of curves y‘=y‘(x), x1: - - - , m, which satisfy 
certain differential equations 

(2) ba(x, (x), ¥(x)) = 0, a=1,---,m<n, 
and certain end-conditions 

(3) Vu(%1, y(%1), X2, y(x2)) = 0, p=i1,---,pS 2nm+2. 


We suppose the usual conditions‘ of continuity and differentiability 


1 Presented to the Society, September 8, 1939. 

2 E. J. McShane, On multipliers for Lagrange problems, American Journal of Math- 
ematics, vol. 61 (1939), pp. 809-819. 

? At the suggestion of Professor M. R. Hestenes. 

4 See, e.g., G. A. Bliss, The problem of Bolza in the calculus of variations, Annals of 
Mathematics, (2), vol. 33 (1932), pp. 261-273. 
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satisfied, as well as the usual conditions on the ranks of the matrices 
of derivatives of ¢. and y,. 

Consider a curve C: y;=yi(x), x1 S42, which is of class C’’ and 
satisfies equations (2) and (3). If X» is a constant and A,(x) a set of 
functions, we define 


F(x, y’; of (x, y’) + y’). 


The statements that C satisfies the multiplier rule, the strength- 
ened Weierstrass condition II and the strengthened Clebsch condi- 
tion III’ shall have their usual meanings. The curve C satisfies the 
transversality condition if there are numbers @, - - - , €¢, such that 


(4) (F — yf Fy,)dx|t + + + = 0 


is an identity in dx, dx2, and 

For normal problems, in which we can suppose \») = 1, the second 
variation of J[y] due to the variations [£:, £, 9(x)] is given by a 
functional ® 


(5) Qlés, n(x), &2, n(x2)] 2w(x, 9, 9’)dx. 


Even for abnormal problems we shall say that the curve C satisfies 
condition IV’ if the functional (5) is positive for every variation 
[£1, 2, n(x) | not identically zero which satisfies the equations of vari- 
ation 


© garni t =0, w= 
where as usual 


My, £2, n2) + yi (x1) + 
+ + yi (x2) + Vuy 


We do not need the specific form of the quadratic function Q or of 2w, 
but for a later paragraph it will be important to note that in (5) the 
integrand f enters only via the combination F(x, y, y’, \), while in the 
coefficients of Q the functions g and y, enter only in the combination 
Nog +e,y, and its derivatives. 


(8) 


5 See, e.g., M. R. Hestenes, Sufficient conditions for the problem of Bolza in the cal- 
culus of variations, Transactions of this Society, vol. 36 (1934), pp. 795, 809, 813. 

* For this computation, see, e.g., G. A. Bliss and M. R. Hestenes, Sufficient con- 
ditions for a problem of Mayer in the calculus of variations, Transactions of this So- 
ciety, vol. 35 (1933); in particular, pp. 311, 312. 
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The sufficiency theorem referred to in the introduction is the fol- 
lowing.” 


THEOREM 1. Jf the curve C is of class C’’ and there are multipliers 
Ao>0, Aa(x) with which it satisfies the multiplier rule, the transversality 
condition and the conditions 11x , 111’ and IV’, then there exists a neigh- 
borhood F of Cin( x, y)-space and a neighborhood N of (x1, y(x1), x2, y(x2)) 
in (x1, 91, X2, Y2)-space with the following properties. If C: y=4(x), 
Sx < lies in F, and the point §(%:), 9(%2)) is in N, and C 
satisfies conditions (2) and (3), then either C is identical with C or 
J[5]>J ly]. 


We now establish a slight generalization of this theorem. 


THEOREM 2. In Theorem 1, the hypothesis \y)>0 can be weakened to 
Ao 20. Moreover, if the other hypotheses of Theorem 1 hold with multi- 
pliers }xo=0, Xa(x), the conclusion can be strengthened by replacing the 
last clause by the words “then C is identical with C.” 


Observe that f, g, Xo never occur in the hypotheses save in the com- 
binations of, Aog and their derivatives. Suppose \» =0. All hypotheses 
continue to be satisfied if we change the problem by replacing Xo 
by 1 and replacing f and g by 0. For this new problem the hypotheses 
of Theorem 1 hold, and clearly the new functional J [y] is identically 
0. Therefore the conclusion of Theorem 1 holds, but the inequality 
J [5|>J[y] is impossible. This establishes Theorem 2. 

Henceforth we assume that the hypotheses of Theorem 2 hold with 
Ao =0. The curve C is then surely not normal; if 


is an aggregate of variations satisfying equation (6), the matrix 


has rank less than p. Let g<p be the greatest rank which the matrix 
(10) can have, and let the set (9) be so chosen that the matrix has 
rank g. By renumbering if necessary, we can bring it about that the 


upper left minor, corresponding to wu, j7=1, - - - , g, is nonsingular. 
Consider now the Mayer problem of minimizing the function 
(11) J* Ly] = y(21), x2, y(%2)) 


7 For a proof, and for reference to other papers on the subject, see M. R. Hestenes, 
A direct sufficiency proof for the problem of Bolza in the calculus of variations, Trans- 
actions of this Society, vol. 42 (1937), pp. 141-154. 
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(where @, - - - , €é, are the constants entering in the transversality con- 
dition) in the class of arcs satisfying the differential equations (2) and 
the end-conditions 


(12) y(x1), %2, y(x2)) =0, q- 


We shall prove that for this problem the curve C is normal and gives 
J* a strong proper relative minimum. That is, the conclusion of Theo- 
rem 1 holds if we replace the reference (3) by (12) and change J to J*. 

The normality of C is evident; for this Mayer problem, the ana- 
logue of matrix (10) is the upper left g-square minor of (10), which 
is nonsingular. 

Equation (11) is the analogue of (1) with 
(13) g* (x1, ¥1, %2, yo) = X2, f*(x, y’) =0. 
Define 
(14) Ae = 1, = 
Then 

F*(x, y, = + ak(x) 9’) 
= ha(x)ba(x, y, = F(x, y, 9’, 

Since the multiplier rule and conditions IIy’ and III’ involve only F 


and its derivatives, they are satisfied for our Mayer problem. 
If we choose 


(16) ec 0, 
we find, recalling (4) and (15), 


(15) 


(F* — yf + lt + + 
2 2 
= (F — yFyy)dx |: + yili + dodg + = 0 
identically in dx,, dxe, dy and dyj2. So the transversality condition is 


satisfied in our Mayer problem. 
From (13), (14) and (16) we have 


(17) Aotg* + dog + 


Now in the light of (17) and (15) the sentence after (8) assures us 
that the sum (5) is identically equal to its analogue for our Mayer 
problem. Hence condition IV’ is satisfied. All the hypotheses of Theo- 
rem 1 have now been verified for our Mayer problem, and it follows 
that C minimizes J* [y], as stated. 
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NEUTRAL ELEMENTS IN GENERAL LATTICES! 
GARRETT BIRKHOFF 


1. Introduction. O. Ore has defined “neutral” elements in modular 
lattices as elements a satisfying an (x u y)=(an x) u (an y) for all 
x, y and dually.? In the case of complemented modular lattices, the 
neutral elements compose the “center” in J. von Neumann’s theories 
of continuous geometries and regular rings—that is, the set of ele- 
ments having unique complements.* 

The purpose of the present note is to extend the notion of neutral 
elements to general lattices. More precisely, call an element a of a 
lattice “neutral” if and only if every triple {a, x, y} generates a dis- 
tributive sublattice. It is proved that the neutral elements of any 
lattice ZL form a distributive sublattice, consisting of the elements 
carried into [J, O] under isomorphisms of L with sublattices of di- 
rect products. Actually, this sublattice is the intersection of the maxi- 
mal distributive sublattices of L. 

Further, complements of neutral elements, when they exist, are 
unique and neutral. The sublattice of complemented neutral elements 
may be called the “center” of a lattice: it consists of those elements 
carried into [I, O] under isomorphisms of L with direct products. 


2. Fundamental definition. We define an element a of a lattice L 
to be “neutral” if and only if every triple {a, x, y} generates a dis- 
tributive sublattice of L. 


Lemna 1. If a is “neutral,” then the dual correspondences xx na 
and x—>x u a are endomorphisms* of L. 


ProoF. By definition, (x u y)n a=(xua)n (yua)and (xn y)n a 
=(xn a) n (yn a), and dually. We note that this condition, which is 
sufficient to guarantee neutrality in the case of modular lattices, does 
not guarantee neutrality in general—see the graph below. 


1 Presented to the Society, September 8, 1939. 

20. Ore, On the foundations of abstract algebra 1, Annals of Mathematics, (2), vol. 
36 (1935), pp. 406-437. For the definitions of lattices and modular lattices (called by 
Ore structures and Dedekind structures), as well as of sublattice, distributive lattice, 
O, I, and so on, compare the author’s Lattices and their applications, this Bulletin, 
vol. 44 (1938), pp. 793-800—or the author’s Lattice Theory, American Mathematical 
Society Colloquium Publications, vol. 25, 1940. 

3 J. von Neumann, Lectures on Continuous Geometries, Princeton, 1935-1936. Cf. 
also R. P. Dilworth, Note on complemented modular lattices, this Bulletin, vol. 45 
(1939), pp. 74-76. 

* We define an endomorphism as a homomorphism of L with itself. 
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Lemma 2. If ais neutral, thenxna=yn aandxua=yua imply 
x=y. 
Proor. By direct computation, using the distributive law twice, 
x= xn(xua) = xn(yua) = (xn y)u(xna) 
= (xn y)u(yna) = yn(xua) = ya(yua) = y. 


Using x in the graph, we see that this condition by itself is also not 
sufficient. However, Lemmas 1-2 together are sufficient to guarantee 
neutrality. 


Fie. 1 


Indeed, consider the correspondence x—[x n a, x ua] from L to 
the product® ST of the sublattice S of elements s<a, with the sub- 
lattice T of elements t2>a. By Lemma 1, it is homomorphic onto a 
sublattice of ST; by Lemma 2, it is one-one; hence it is isomorphic. 
Moreover x—[a, a] = [J, O], since a is the I of S and the O of T. 

But conversely, (7, O] is obviously “neutral” in ST, since each 
component is. Hence it is neutral in every sublattice of ST, including 
L, and we conclude® that the following holds: 


THEOREM 1. An element of a lattice is neutral if and only if it is 
carried into [I, O] under an isomorphism of the lattice with a sublattice 
of a direct product. 


3. Neutral elements a sublattice. Just as in the case of modular 
lattices, we have the following theorem: 


THEOREM 2. The neutral elements of any lattice form a distributive 
sublattice. 


Proor. Let a and 6 be neutral. Then since the product of two 


5 By the “product” ST, is meant the system of couples [s, tl, se5S,teT, where 
[s, t] [s’, t’']=[sms’, tnt’] and dually. 
6 N.B., we do not assume that L itself has an O or an I. 
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endomorphisms of L is an endomorphism, it follows that the cor- 
respondence 


x— [xnanb, (xua)nb, (xna)ub, xuaub] 


defines a homomorphism of L onto a sublattice of a product STUV 
of four sublattices of itself. Moreover by Lemma 2, (x n a) and (x u a) 
and therefore x, are determined uniquely by their images under the 
homomorphism. Hence the homomorphism is an isomorphism. But 
au b goes into [J, I, I, O] under this—and hence into [J, O] of 
(STU) V. Thus au is neutral; dually, a is neutral, which was 
to be proved. 


4. Intersection of maximal distributive sublattices. The set of neu- 
tral elements of a lattice is also characterized in another way by the 
next theorem: 


THEOREM 3. The set of neutral elements of a lattice L is the intersec- 
tion of its maximal distributive sublattices. 


Proor. First, if a is not neutral, then some triple {a, x, y} is not 
distributive. Hence no maximal distributive sublattice obtained by 
enlarging the distributive sublattice generated by {x, y} can con- 
tain a. Consequently, the intersection of the maximal distributive 
sublattices of ZL contains no non-neutral elements. 

Conversely, if a is neutral, and S is a distributive sublattice of L, 
consider the sublattice generated by {a, S}. The endomorphisms 
x—x Na and x—x ua carry it into sublattices generated by a dis- 
tributive sublattice and J or O. But such sublattices are always dis- 
tributive—hence so is the sublattice generated by {a, S}, since it is 
a sublattice of a product of distributive lattices. Thus every maximal 
distributive sublattice contains a, and the intersection of the maxi- 
mal distributive sublattices contains every neutral element (as well 
as no non-neutral elements). 


5. Center of a lattice. When one comes to complements of neutral 
elements, one finds that the following statement holds: 


THEOREM 4. Complements of neutral elements, when they exist, are 
unique and neutral. 


Proor. Using Theorem 1, we see that [Z, O] 9 [x, y] is [O, O] if 
and only if x =O, while [Z, O] u [x, y] is [J, I] if and only if y=TJ. 
Hence [J, O] has no complement except [O, I] in the sublattice of ST 
isomorphic with L, proving uniqueness. Moreover [O, I] is itself neu- 
tral, completing the proof. 
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Coro.iary 1. The neutral elements of a complemented lattice form a 
Boolean algebra. 


We define the “center” of a lattice as the set of its complemented 
neutral elements. 


THEOREM 5. The center of any lattice L is a complemented distribu- 
tive sublattice—and hence a Boolean algebra. 


ProoF. If a and b are neutral elements of L, with (neutral) com- 
plements a’ and 0’, then 


(amb) a (a’u = (anbna’)u (anbnd’) =Ou0 =O, 
(anb)u (a’u = (aud ud’) a(bud ud’) =InI=T, 


and so an bd is complemented. Dually, a u 6 is complemented, com- 
pleting the proof. 
We can now specialize Theorem 1 by proving the following: 


THEOREM 6. An element ts in the center of a lattice L if and only if 
it is carried into [I, O| under an isomorphism of L with a direct product. 


Proor. By Theorem 1, such an element is neutral, and it has the 
complement [O, I]. Conversely, suppose a and a’ are complementary 
neutral elements of L. Then for all x, 


x= xnIT=xn(aua’) = 


Hence the correspondence x—>[x na, xn a’] is one-one between L 
and the couples [u,v] with u<a, v<a’; the inverse correspondence 
is [u, v]—u vu v. But it obviously preserves inclusion; hence it is an 
isomorphism. Finally, it carries L into the product ST of the lattice S 
of elements s <a, with the lattice T of elements ¢ <a’, while it carries a 
into [ana, ana’|=[I, O] in ST. 
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ON TOPOLOGICAL COMPLETENESS' 
L. W. COHEN 


Recently A. Weil? defined a uniform space as a set of points p such 
that for each a in a set A there is defined a set U,(p) ¢ S, the class 
of sets U.(p) satisfying the conditions: 

II. U.(p) = (p). 
II,. To each a, B 2 A there is a y=7(a, B) e A such that U,(p) 
c U.(p) Us(p). 

III,. Toeach ae A there is a B(a) A such that if p’, p’’ (Q), 
then p’’ e U,(p’). 

For the uniform space S, Weil introduced the concept of Cauchy 
family {Mg} of sets. Such a family is defined by the conditions that 
the intersection of any finite number of sets of the family is non- 
empty and that to each a e A there is a p, e S and a B(a) such that 
Mga) ¢ Ua(pa). Weil gives a theory of completeness in these terms. 

The writer has considered* a space S of points p and neighborhoods 
U.(p) where a@ is an element of a set A such that: 

I. [].U.(p) =(0). 

II. To each a and Be A and pe S there is a y=y(a, B; p) such 
that U,(p) ¢ Ua(p) Usp). 

III. To each ae A and pe S there are X(a), 5(p, a) e A such that, 
if Uxcay(P) £0, then Us.y,«)(g) Ua(P). 

The uniformity conditions here are lighter than those in II, and 
III,. A Cauchy sequence p, ¢ S was defined by the condition that for 
every ae A, n, and p, S exist such that p,2¢ U.(pa) for S is 
complete if every Cauchy sequence has a limit. It was shown that 
there is a complete space S* which contains a homeomorphic image 
of S such that the image of a Cauchy sequence in S is aconvergent 
sequence in S*. 

It is the object of this paper to show that Weil’s space is a special 
case of the space S;,11,1 and that the notion of Cauchy family in this 
space leads to the same theory of completeness as that previously de- 
veloped. 


THEOREM 1. Jf S satisfies III, and B?(a)=B(B(a)), then from 
Us (a)(q) (Pp) £0 follows ¢ Ua(p). 


1 Presented to the Society, December 27, 1939. 

2 A. Weil, Sur les Espaces & Structure Uniforme, Paris, 1938. 

?L. W. Cohen, On imbedding a space in a complete space, Duke Mathematical 
Journal, vol. 5 (1939), pp. 174-183. Also Duke Mathematical Journal, vol. 3 (1937), 
pp. 610-617, where the notion of topological completeness is introduced. 
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Proor. Let s ¢ Us (a)(q) Usa) (p). Then from s, g’ Us we 
have q’ ¢ Ug,a)(s). Therefore ¢ Similarly Ug a)(p) 
Ugia)(s). Now from p Ugiay(s) and Us a)(g) ¢ we have 
Upa)(q) ¢ 


Coro.iary. S satisfies 1114, then S satisfies 


Proor. For any p ¢« S and a e A we need only take X(a) = 5(p, a) 
=6(8(a)). The result is stronger than III since 5(p, a) is independent 
of p. 

From now on a space S is one satisfying I, II, III. A family of sets 
{ Mg} is a Cauchy family if the intersection of any finite number of 
Mz is non-empty and if for any aeA there is a B(a) such that 
Mga) ¢ Ua(pa) for some pa e S. We will say that S is W-complete if, 
for every Cauchy family {Mg}, [[sM@s~0, where Mz is the closure 
of Mg. We shall always use the notations (a), 5(p, a) in the sense 
of III. 


THEOREM 2. S is W-complete if and only if every Cauchy family 
{ Ua(pa)} consisting of one U2(pa) for each a e A has the property that 
#0. 


Proor. Assume S is W-complete. If { Ua(pa)} is a Cauchy family, 
then #0. Since ¢ Ua(p) follows from III, the 
condition [].Us(Pxa)) #0 is necessary. Assume now that the condi- 
tion is satisfied and let {M,} be a Cauchy family in S. To each 
aeA there are p, e Sand Mga) such that Mga) ¢ Ua(pa). It is clear 
that {U.(p.)} is a Cauchy family. Hence there is p e []«Ue(Pr.a)- 
For any yeA, consider the A(y) and a=46(p, y) of III. Since 
Ual Pray) Ual Pray) Uy(p) and 


© © Val Prcay) 
Hence 


for all 8, y. Thus p and S is W-complete. 

The space S* referred to above is defined as follows. A family 
{ Ua(pa)}, one for each a e A, such that for any (a1,---, @) CA, 
Ua,(Pa;) Ua,(Pa,) #0 is denoted by II. We write II’~II”’ 
if for every a e A there is a set (a1, -- - , @,) CA such that for some 


pate S 
II + I € Vale): 
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The relation II’~II”’ classifies all II into mutually exclusive classes 
C(I) which are the points P of S*. The neighborhoods U%,...«,(P) 
are made up of all Q=C(II®) S* where II? = { Ua(ga) } is such that 
for some Qn; P), t=1,---, m=m(a,---, Qn; P), 


II II Ua,(Pa;)- 
i=1 i= 
The space S* (for which A* is the set of all finite subsets of A) satis- 
fies I, II and 
III*. For each (a, ---, @,)¢A and for Pe S* there exist 
+>, Qn; P), #=1,---, m=m(m,---, Gn; P), such that 
The mapping f(p) = C(II”), where II? = { U.(p) } , is the homeomor- 
phism on S to a subset of S* referred to above. 


THEOREM 3. S* is W-complete. 


Proor. We first show that if { U%,....,(P"°°'™)}, where each 
(a;, - - - , @,) © A occurs just once, is a Cauchy family in S*, then 


Il #0. 
Consider the U.*(P*) for alla eA. Now P*=C(II"*), = { U,(p%)}, 
U,(p$) being a neighborhood in S for each (a, y) ¢ A. Since, for each 
(a1, Qn) cA, 0, we have 


II 0. 

i=1 
Hence the family II = { U.(p2) } defines a P=C(II) e S*. For any 
(1, +Y¥m) CA 


II U3AP) = 

i=l 
Both U%,(P) aud U*,(P%) are the set of all Q = C(II2), 11? = { Us(qa) } 
such that for some By =B,(7:),j= 1,---,ki, 
Hence 


Ut (P) = (P) = II (Pr), 


Thus from the fact that { U%,...2,(P™°° ‘an)} is a Cauchy family we 
have for any sets (a1, - ,@n), CA 


— 
— 
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It follows that 
(1) Pe 


(ay+++a,)CA 


Now suppose that {Ms} is a Cauchy family of sets in S*. For 
every (V1,°°-°, Yn) CA there is a B(j1, --- , Yn) such that for some 
S*, MB © Thus 


i=1 i=1 


for every finite set (yi,---, and { isa 
Cauchy family. Let P ¢ S* satisfy (1). For any (a, --- ,a@,) ¢ A and 
P e S* let Ym) CA be the set of - - - , @n; P) satis- 
fying III*. From III* and (1) 


Us,..-a,(P). 
Since { M#} is aCauchy family, for any Band (ay, - - - ,a,) ¢ A we have 


¥ 0. 
Hence P e[][sM# and S* is W-complete. 


THEOREM 4. If { Mg} is a Cauchy family in S and f(S) ¢ S* is the 
homeomorphism defined above, then there is a P ¢ S* such that 


1. [[sf(Ms) =(P), 
2. for any (a1, ---,Q@n) CA there are - - - , Bm such that 


U%,...»,(P). 


i=l 


ProorF. f(Ms,) is the class of P = C(II”) for all » « Mg where for any 
(a1, CA and there are B(A(a;)) and such that 


Ma ca;)) © (Prca;) 


i=1 i=1 
If ge] then for 6;=8(pr(@;), we have 
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qe II U;,(9) II Val 


t=1 t=1 


the family { U.(Pria)) } =II defines a P = C(I) e S* such that 


Hence for such (a1,---, @,)CA there are and 
- - - , m, such that 


i=l 


Thus the family of all finite products {f(Mz,) -f(Mz,)} is a Cauchy 
family in S* since f(Msg,) - - - f(Ms,) =f(Mz, - - - #0. We have 


i=1 


for all 8 and all (a, - -- ,a@,) ¢A. In other words 
Pel f(Ms). 


But the intersection of all f{( Msg) is P, for if P’ is any other point there 
are sets (ay, ---,@n), (a1, °° such that 


= 0 
and 


P’ eS* — — — Tf). 
B 


i=1 


We conclude with the remark that if S is W-complete, S is com- 
plete. Suppose p, is a Cauchy sequence in S. Let Mn=(pn, Pau, + *)- 
Then the intersection of any finite set of M, is non-empty and for 
any a¢A there is an m(a) such that My,2) ¢ Ua(pa) for some pa e S. 
Thus {M,} is a Cauchy family. S being W-complete, there is a 
pe]],M,. Now for p, any ae A, \(a) and 6 = 6(p, a), we have 


Ms) ¢ Us( ps), 0 Uy(a)(p)M acs) Urcay(p) Usps), 
Mas) Us( ps) U.(p), Pn U.(p), n n(5(p, lim Pn =p. 
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